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Preface 


The  study  of  aircraft  flight  control  is  a  very  broad  field.  In 
selecting  a  thesis  topic,  I  was  interested  in  choosing  a  project  that 
would  permit  me  to  experience  as  many  aspects  of  this  field  as  possible 
This  necessitated  that  the  project  be  small  enough  in  scope  that  it 
afforded  me  the  opportunity  to  engage  in  analysis  as  well  as  design. 

The  NEXTRPV  program  was  well  suited  to  meet  this  goal.  Whereas  this 
project  is  not  necessarily  glamorous,  it  gave  me  the  chance  to  exper¬ 
ience  "basic"  engineering  which  I  believe  is  an  important  part  of  any 
education. 

This  thesis  originally  began  as  strictly  a  flight  control  system 
analysis  and  design  for  an  existing  remotely  piloted  vehicle  (RPV).  As 
the  project  progressed,  it  became  evident  that  the  RPV  was  fairly  stabl 
and  would  not  need  as  much  flight  control  augmentation  as  initially 
anticipated.  Therefore,  this  thesis  was  extended  to  include  an  invest! 
gation  of  a  practical  application  of  this  RPV. 

Many  people  contributed  to  this  thesis.  To  each  of  them,  I  say 
thank  you.  Special  thanks  goes  to  my  thesis  advisor,  Qaptain  James  T. 
Silverthorn,  for  his  numerous  suggestions  and  neverending  attention. 

My  gratitude  is  extended  to  Professors  Robert  A.  Calico  and  John  J. 
D'Azzo  whose  assistance  was  invaluable.  My  appreciation  is  also  ex¬ 
pressed  to  Mrs  Anna  Lloyd  for  preparing  this  manuscript.  Lastly,  I 
would  like  to  thank  my  wife,  Susan,  whose  constant  support,  encourage¬ 
ment,  and  understanding  made  the  completion  of  this  project  possible. 

Brian  L.  Jones 
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Abstract 


~  ^Tbis  study  investigated  the  stability  and  control  of  a  remotely 
piloted  vehicle  (RPV)  with  a  thrust  vectoring  unit  attached.  All 
geometric  and  aerodynamic  data  was  generated  and  used  to  analyze 
the  RPV.  Specific  handling  qualities  were  developed  and  compared 
with  the  RPV  characteristics.  This  comparison  indicated  that  the  RPV 
was  too  oscillatory  in  both  the  phugoid  and  dutch  roll  modes.  Also, 
the  RPV  displayed  a  dominant  spiral  mode.  A  flight  control  system 
was  synthesized  to  eliminate  those  traits.  Evaluation  of  this  flight 
control  system  was  conducted  through  the  use  of  three  different 
sensitivity  studies  and  a  nonlinear  simulation.  In  addition,  a  model 
matching  application  was  examined  for  this  RPV,  Model  matching  entails 
using  design  procedures  to  synthesize  an  expanded  flight  control  system 
so  that  the  RPV  has  dynamic  characteristics  similar  to  the  F-15. 
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FLIGHT  COriTilOL  SYSTEM  ANALYSIS  AND  DESIGN  FOR 
A  REMOTELY  RILOIED  VEHICLE  WITH  THRUST  VECTORING  UNIT 

I.  Inf.rodiiction 


Background 

In  aircraft  design,  the  aeronautical  engineer  is  faced  with  many 
design  objectives.  For  military  aircraft,  these  objectives  would  in¬ 
clude  such  traits  as:  low  observability,  high  lethality  and  survivabil¬ 
ity,  STOL  capability,  and  lastly,  air-to-air  and  air-to-ground  capabil¬ 
ity.  Several  schools  of  thought  exist  on  how  to  satisfy  thci^e  objectives. 
One  alternative  is  the  use  of  a  thrust  vectoring  unit  (TVU).  With  such 
a  unit  employed,  the  tail  section  could  be  removed  without  any  reduction 
in  aircraft  control.  This  would  reduce  the  infrared  signature  and  hence 
lower  an  aircraft's  observability.  Also,  thrust  vectoring  ificreases  an 
aircraft's  maneuverability,  which,  in  turn,  increases  its  lethality  and 
survivability  while  enhancing  its  air-to-air  and  air-to-ground  cooabil- 
ity.  In  addition,  thrust  vectoring  has  been  shown  to  improve  an  air¬ 
craft's  STOL  capability. 

With  a  continual  reduction  in  funding,  the  Air  Force  needs  to  find 
inexparisi VG  ways  in  which  to  develop  and  test  such  theories  as  thrust 
vectoring.  One  answer  to  this  pr'^blem  is  using  reicotely  piloted  ve¬ 
hicles  (RPVs)  as  flight  test  vehicles.  RPVs  can  be  operated  at  a  small 
fraction  of  the  cost  of  a  test  aircraft,  resulting  in  substantial 
savings  in  research  dollars. 

The  two  above  mentioned  concepts  were  combined  to  foi'in  the  NEXIRPV 
program  under  the  direction  of  the  Control  Systems  Dove! upniont  Branch 


(AFWAL/FIGL ) ,  Flight  Dynamics  Laboratory,  Air  Force  Wright  Aeronautical 
Laboratories.  Their  stated  objective  was  to 

demonstrate  the  use  of  practical  thrust  vectoring  and  digital 
flight  control  systems  in  a  low  cost  flight  research  vehicle 
oriented  toward  a  more  stealthy  design  with  an  expanded  maneuver¬ 
ing  envelope. 

The  approach  the  laboratory  intended  to  pursue  was  to  first  integrate 
the  TVU  into  the  RPV  through  an  analog  control  system.  Flight  control 
authority  would  gradually  be  transferred  to  the  TVU  over  a  series  of 
test  flights.  Following  the  successful  completion  of  the  TVU  full 
authority  test  flight,  the  entire  procedure  v;ould  be  repeated  using  a 
digital  control  system.  At  the  conclusion  of  this  phase,  the  tail 
section  would  be  removed  for  a  final  performance  evaluation. 

The  RPV  chosen  for  this  program  suited  it  well  for  several  reasons. 
Since  this  RPV  was  obtained  from  a  previous  program,  it  had  already  been 
flight  proven.  Thus,  all  considerations  other  than  flight  control  system 
(FCS)  Oesign  were  eli.ninated  from  this  project.  Also,  the  twin  boom/ 
vertical  tail  assembly  required  little  modification  for  the  TVU  addition 
(see  Fig  ").  Lastly,  since  the  RPV  v/as  already  constructed,  the  associ¬ 
ated  develcpiiient  costs  were  saved. 

Problfm  Statement <_^op£,_and  Assiiruptions 

The  purpose  of  this  study  v;as  to  determine  the  stability  and  con¬ 
trol  characteristics,  develop  specific  handling  qualities,  synthesize  a 
FCS  to  satisfy  these  handling  qualities,  and  investigate  a  practical 
application  of  this  RPV. 

This  study  was  restricted  to  the  TVU  zero  authority  case.  In  other 
words,  the  TVU  was  attached  to  the  RPV,  but  was  given  no  control 
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authority.  The  only  exception  appeared  in  the  basic  analysis.  In  this 
phase,  center  of  gravity  location  and  TVU  attachment  v/ere  used  as 
parameters.  X-axis  center  of  gravity  (xcg)  assumed  three  locations: 

10,  25,  and  40  percent  mean  aerodynamic  chord  (mac).  The  TVU  attached 
and  unattached  cases  were  examined  for  each  xcg  location.  The  practical 
application  was  restricted  to  the  longitudinal  case.  As  shown  in  Chapter 
V,  the  lateral -directional  case  would  simply  be  a  repetition  of  the 
same  theory. 

Initially,  several  simplifying  assumptions  were  imposed.  Perturba¬ 
tions  about  straight  and  level  flight  were  assumed.  This  allowed  the 
linearized  aircraft  equations  of  motion  to  be  decoupled  into  longitudinal 
and  lateral-directional  sets.  For  the  FCS  design,  only  pitch,  roll,  and 
yaw  rates  were  acceptable  feedbacks.  It  was  further  assumed  that  a 
processing  unit  would  be  accessible  should  compensation  be  required. 
Lastly,  full  state  feedback  was  assumed  available  for  the  practical 
application. 

Approach  and  Presentation 

This  study  was  separated  into  four  major  phases.  Phase  one  con¬ 
sisted  of  obtaining  all  necessary  geometric  and  aerodynamic  data  for 
this  RPV ,  follov.'ed  by  the  generation  of  six  sets  of  stability  derivatives 
(TVU  attached  and  unattached  cases  for  three  xcg  locations).  These  sets 
made  up  the  mathematical  models  used  in  this  study.  Chapter  II  of 
this  report  presents  the  results  ot  phase  one. 

The  basic  analysis  of  the  RPV  constituted  phase  two.  This  analysis 
began  with  the  development  of  specific  RPV  handling  qualities.  The 
aircraft  equations  of  motion  were  derived  and  a  modal  analysis  of  the 
six  models  generated  in  phase  one  v/as  conducted.  Phase  tv/o  concluded 
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with  a  comparison  of  the  modal  analysis  with  the  RPV  handling  qualities. 
The  report  on  phase  two  is  contained  in  Chapter  III. 

Phase  three  involved  the  design  and  evaluation  of  a  PCS  to  satisfy 
all  the  handling  qualities.  Sensitivity  studies  and  a  nonlinear  simu¬ 
lation  were  used  to  evaluate  the  PCS.  Phase  three  is  documented  in 
Chapter  IV. 

Lastly,  phase  four  entailed  an  investigation  into  a  practical  appli 
cation  (model  matching)  for  this  RPV.  An  expanded  PCS  was  synthesized 
to  make  this  RPV  have  dynamic  characteristics  similar  to  an  P-15.  Three 
different  design  procedures  were  examined  and  are  described  in  Chapter  V 

The  final  chapter  of  this  report,  Chapter  VI,  presents  conclusions 
and  recommendations. 
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1 1 .  RP^  Generation 

Aircraft  stability  and  control  characteristics  are  usually  examined 
based  upon  a  linearized  set  of  differential  equations  of  motion.  The 
development  of  these  equations  (see  Chapter  III)  contains  numerous  param¬ 
eters  which  fall  into  three  categories:  physical  quantities,  flight  con¬ 
dition  defining  quantities,  and  stability  derivatives.  The  purpose  of 
this  chapter  is  to  generate  these  parameters  for  this  particular  RPV. 

Method 

The  physical  quantities  (mass,  moment  of  inertia,  etc.)  and  flight 
condition  defining  terms  (airspeed,  dynamic  pressure,  etc.)  v/ere  easily 
determined  from  the  size  of  the  RPV  and  its  expected  flight  regime. 

Thus,  the  major  crux  of  the  model  generation  becam.e  calculating  the 
stability  derivatives.  Applicable  publications  were  consulted  in  order 
to  determine  the  most  appropriate  method  of  calculating  these  deriva¬ 
tives.  Most  publications  reviewed  listed  The  USAF  Stability  and  Control 
DATCOM  (Ref  8)  as  a  reference.  These  manuals  contain  empirical  formulas 
and  methods  for  determining  stability  and  control  derivatives,  aircraft 
moments  of  inertia,  etc.  A  computerized  version  of  these  methods  appears 
Hie  USAF  Stabil i ty  and  Control  Digi tal  Datcom  (Ref  18) .  Due  to  the 
widespread  acceptance  of  Datcom  methods  and  the  versatility  of  a  com¬ 
puter  program,  it  v/as  decided  that  Digital  Datcom  would  be  the  best 
source  of  stability  derivatives. 

Input  Data  fqji  Digi  tal  Datcom 

Digital  Datcom  requires  three  types  of  input  data:  aircraft 
geometry,  airfoil  section  characteristics,  and  flight  condition. 


Aircraft  geometry  was  obtained  by  directly  measuring  the  RPV's  physical 
dimensions.  A  circular  cross-section  approximation  was  employed  for 
the  fuselage  and  booms  (see  Fig  1).  Table  I  is  a  summary  of  all 


Figure  1.  Circular  Cross-Section  Approximation 
for  tfie  Fuselage  and  Boom 


measurements. 

Two  major  modifications  were  planned  for  the  RPV  and  were  incor¬ 
porated  into  the  mathemacical  model.  The  first  was  a  two  foot  nose 
extension  and  ll:e  second  was  the  addition  of  the  TVU.  Both  modifica¬ 
tions  were  modelled  as  part  of  the  fuselage. 

Next,  airfoil  section  characteristics  were  required  for  each  air¬ 
foil.  The  wing  was  found  to  have  a  va-ying  airfoil  section,  a  NACA 
Z3012  at  the  root  and  a  NACA  44i2  at  the  tip.  The  horizontal  tail  was 
a  NACA  0009  while  the  vertical  tails  were  NACA  0010  airfoil  sections. 
Since  NACA  airfo'’!  sections  w^re  found  (Ref  1),  Digital  Datcom's  in¬ 
ternal  aerodynamics  package  v/as  utilized  for  the  generation  of  section 
characteristics.  An  appreximation  of  the  wing's  varying  airfoil  section 
was  formulated  by  computing  two  sets  of  stability  derivatives,  one  set 
for  each  of  the  root  and  tip  airfoil  sections,  and  linearly  averaging 
the  results.  This  approximation  should  be  accurate  since  the  differ¬ 
ence  between  the  individual  stability  derivatives  was  small. 
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Table  I.  Physical  Din’cnsioris  of  the  RPV 
a.  Wing,  Horizontal  Tail,  and  Vertical  Tail  Dime 


W. 

H.T. 

Apex,  ft: 

x-axis 

4.25 

11.39 

z-axis 

0.60 

2.10 

Chord,  ft: 

root 

2.52 

1 .29 

tip 

2.00 

1 .29 

Incidence  Angle,  deg: 

0.95 

0 

Span,  ft: 

exposed 

13.00 

4.88 

theoretical 

13.71 

5.00 

Sweptback  Angle,  deg: 
leading  edge: 

inboard  panel 

3.50 

0 

outboard  panel 

3.50 

0 

trailing  edge: 

inboard  panel 

0 

0 

outboard  panel 

-1 .40 

0 

*  See  Ref  18  for  definitions  of  terms. 
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Table  I.  (continued) 


b.  Fuselarie  Dimensions 

j 

‘  x-station,  ft 

perimeter,  ft 

9 

area,  ft 

radi us , 

2.00 

2.42 

2.83 

3.25 

3.67 

1  3.83 

i  4.25 

i  4.75 

5.25 

5.75 

6.25 

6.75 

3.33 

3.42 

3.58 

3.71 

3.79 

3.83 

3.75 

3.75 

3.75 

3.75 

3.75 

3.75 

0.88 

0.93 

1.02 

1  .09 

1.14 

1.17 

1.12 

1.12 

1.12 

1.12 

1.12 

1.12 

0.53 

0.54 

0.57 

0.59 

0.60 

0.61 

0.60 

0.60 

0.60 

0.60 

0.60 

0.60 

c .  Boom 

Dimensions 

1 

I  x-station,  ft 

perimeter,  ft 

area,  ft^ 

radius. 

ft 

5.00 

0.67 

0.035 

o.n 

5.42 

1 .58 

0.20 

0.25 

5.83 

1 .83 

0.27 

0.29 

6.33 

1.92 

0.29 

0.30 

6.75 

1.67 

0.22 

0.26 

7.17 

1 .62 

0.21 

0.26 

7.58 

1.58 

0.20 

0.25 

8.00 

1.50 

0.18 

0.24 

1  8.42 

1 .42 

0.16 

0.22 

!  8.83 

1 .35 

0.15 

0.22 

9.25 

1.29 

0.13 

0.21 

1 

9.67 

1.21 

0.12 

0.19 

1 

1  10.03 

1.12 

0.10 

0.18 

10.50 

1  .00 

0.080 

0.16 

10.92 

0.92 

0.067 

0.15 

11.25 

0.83 

0.055 

0.13 
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Flight  condition  defining  paracietors  v;as  tr  n  loot  oet  of  input  data 
needed.  After  consideration  of  the  engine  capabilities  and  the  flight 
profiles  expected,  the  airspeed  regime  was  determined  to  range  from 
Nach  0.05  to  Mach  0.20.  Similarly,  the  angle  of  attack  range  was 
formulated  as  -4  deg  to  15  deg.  Sea  level  conditions  were  assumed 
for  all  quantities  requiring  atmospheric  conditions. 

Formul ation  U]e_  R^V  Representations 

The  choice  of  Digital  Datcoin  specified  the  basic  RPV  representation 
as  a  straight  tapered  wing  and  one  vertical  tail  of  ircr'-^ased  area.  It 
was  decided  to  employ  a  scale  factor  of  1.75  in  increasing  the  vertical 
tail  area.  This  choice  was  somewhat  arbitrary  since  no  documentation 
could  be  found  on  this  type  of  approxination.  Geometric  relationships 
(sweepback  angles)  and  aspect  ratio  wore  preserved  in  the  scaling  up  of 
the  vertical  tail. 

In  addition,  the  cross-sectional  areas  (in  the  body  namelist)  were 
increased  to  include  the  effects  of  the  booms.  Because  of  the  blanking 
effect  of  the  TVU,  only  the  area  of  one  boom  was  added  to  the  fuselage. 
When  the  TVU  was  absent  from  the  RPV  configuration,  tfie  areas  of  both 
booms  were  added  to  the  fuselage.  All  other  data  util  iced  in  the 
representations  was  taken  from  Table  1. 

Two  par-ameters  were  chosen  for  investigation:  xcg  location  and 
TVU  attachment.  The  xcg  location  assumed  three  values  (10,  25,  and 
40  percent  mac).  Goth  the  TVU  attached  and  unattached  cases  were 
examined  for  each  xcg  location.  Thus,  a  total  of  six  cases  were 
investigated: 
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Case  1  -  TVU  attached,  xcg  -  10?^  mac 

Case  2  -  TVU  attached,  xcg  =  25«  mac 

Case  3  -  TVU  attached,  xcg  =  mac 

Case  4  -  TVU  unattached,  xcg  =  10'’^  mac 

Case  5  -  TVU  unattached,  xcg  =  mac 

Case  6  -  TVU  unattached,  xcg  =  40%  mac 

Hence,  the  effect  of  the  TVU  and  the  xcg  location  on  the  RPV  character¬ 
istics  could  be  detennined. 

With  the  input  data  finalized,  the  RPV  representations  were  com¬ 
plete.  Figure  2  contains  the  input  data  cords  for  Digital  Datcom. 

The  four  underlined  statements  had  to  be  changed  with  each  case.  A 
visual  interpretation  of  the  RPV  with  the  TVU  attached  and  unattached 
can  be  found  in  Figs  3  and  4,  respectively. 

Data  Reduction 

Since  the  entire  flight  regime  specified  16  Mach  nun-bers  and  20 
angles  of  attack  per  Mach  number.  Digital  Datccr.i  generated  320  sets  of 
stability  derivatives  for  each  case  (i.e.,  a  data  bank).  Unfortunately, 
the  aircraft  equations  of  motion  could  accept  only  one  set  of  stability 
derivatives  at  a  time.  Thus,  a  nominal  flight  condition  was  needed. 
Three  techniques  were  examined  in  the  selection  of  this  flight  con¬ 
dition:  lift  averaging,  best  endurance,  and  best  range. 

The  lift  averaging  technique  assumed  straight  and  level  flight. 

Thus , 

W  =  L 

-  I  pV^SCl 
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3.  Ccnipririscn  of  the  Actual  RPV  and  Its  Di 
Datcoin  Representation  with  the  TVU  Attached 


Figure  4.  Ccti'p.irison  of  the  Actu.il  RPV  and  Its  Digital 
Datcoiii  Representation  with  the  TVU  Unattached 
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Substituting  and  solving  for  Cl  yields 


CL 


Spc2m2 


Assuming  the  n.aximum  weight  (250  lb)  and  substituting  the  appropriate 
values  gives 


or 


Cl 


5.45  X  10-3 
m2 


•30  J<J0li 

Cl1/2 


(1) 

(2) 


From  Eq  (1),  Cl  can  be  calculated  over  the  range  of  reasonable  Mach 
numbers  (Table  II).  A  linear  average  of  lift  coefficient  was 
determined  as  Cl  =  0.59  which  represented  the  average  amount  of  lift 
required  over  the  entire  flight  regime.  Substituting  this  value  into 
Eq  (2)  produced  an  associated  Mach  number  of  Mach  0.10.  Digital  Datcoin 
specifies  its  flight  condition  in  terms  of  angle  of  attack  rather  tiian 
the  lift  coefficient.  Therefore,  the  data  bank  v/as  interpolated  (using 
Mach  0.10  and  Cl  =  0.59)  to  produce  an  associated  angle  of  attack  of 
a  =  2  deg. 

As  alternate  methods  for  selecting  a  nominal  flight  condition, 
best  endurance  and  best  range  flight  profiles  v/ere  considered  (Ref  7: 
Sec  III,  33,  37).  After  considering  that  the  RPV  would  be  flovai  oy 
sight,  it  became  apparent  that  neither  of  those  two  flight  profiles 
could  be  maintained.  Tnerofore,  the  results  of  the  lift  averagino 
technique  (Mach  0.10,  -  2  deg)  wore  used  as  the  nominal  fligiit 

condi tion. 


14 


Finally,  the  trim  condition  needed  to  be  verified  to  insure  that 
this  was  an  appropriate  choice  of  flight  condition.  During  straight  and 
level  flight,  the  pitching  moment  equation  becomes 

Cm  ^  Ciuq  +  a  +  6e 

In  a  trimmed  aircraft.  Cm  =  0.  Using  this  value  and  solving  for 
produces 

^Gtnin  -  - — 

Cfiioe 


Substituting  the  Case  2  values 


■istr’ini 


-0.06  +  (0.023)(2.0) 
■  -07034 


=  0.33  deg 


Thus,  it  was  decided  that  the  RPV  was  in  a  trimmed  flight  condition. 
Table  III  contains  a  complete  summary  of  the  stability  derivatives  at 
the  nominal  flight  condition. 


rt^d  Calculations  of  liomiinin^  Terms 

Three  types  of  parameters  were  still  needed  to  complete  the  model 
of  the  RPV:  the  remaining  stability  derivatives  (e.g.  Cx,  J ,  the  tontrol 
derivatives  (e.g.  Cit;.^),  and  the  physical  quantities  (e.g.  ly) .  The 
complete  derivation  of  these  parameters  is  given  in  Appendix  A.  In 
addition,  these  parameters  are  siimmarized  in  Table  III. 


Summary 


Table  III  contains  a  summary  of  the  mathematical  models  used  in 
this  study.  Five  approximations  were  included  in  the  generation  of 
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nd  ccr.trcl  derivatives  have  units  of  per  rad. 


these  models; 


1.  Circular  cross-sections  were  used  for  the  fuselage. 

2.  A  straight  tapered  wing  was  substituted  for  the  cranked  wing. 

3.  One  vertical  tail  of  increased  area  (75%  greater)  was  sub¬ 
stituted  for  the  twin  vertical  tails. 

4.  The  effects  of  the  booms  were  incorporated  into  the  fuselage 
representation. 

5.  Average  stability  derivatives  were  computed  using  the  root  wing 
section  and  tip  wing  section  cases. 
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III.  Analysis  of  the  RPV 


Before  any  roasonable  PCS  synthesis  could  bo  undertaken,  flight 
specifications  had  to  be  developed  and  the  unaugmented  system 
analyzed.  Thus,  the  purpose  of  this  chapter  is  to  present  specific 
RPV  handling  qualities,  followed  by  the  above  laentioned  analysis. 

This  analysis  was  divided  into  four  parts:  derivation  of  the  aircraft 
equations  of  nioticn,  examination  of  the  modal  ciiai'actoristics ,  genera¬ 
tion  of  transfer  functions  and  time  responses,  and  development  of  con¬ 
clusions. 

Handl i ng  Qual  i  ties 

A  literature  review  resulted  in  only  one  publication  (Ref  14)  which 
specifically  addressed  RPV  handling  qualities;  unfortunati^ly ,  sections 
of  this  publication  are  incomplete.  However,  since  all  publications 
reviewed  relied  heavily  upon  the  criteria  set  forth  in  various  Military 
Specifications,  it  was  decided  that  a  comparison  study  of  the  following 
four  specifications  would  provide  the  most  appj'opriate  sout'co  of  liandling 
qualities: 

1.  RPV  Flying  Qualities  Design  Criteria,  A'TDL -TR-yC-l 25  (Ref  14). 

2.  Flying  Qualities  of  Piloted  Aircraft,  MlL-F-o/BbB  (Ref  4). 

3.  Flying  Qualities  of  Piloted  V/STOL  Aircraft,  MIL-F-83200  (Ref  5). 

4.  Flight  Control  Systems  -  Design,  Instal la tion ,  and  Test  of 
Piloted  Vehicles,  niL-F-y4900  (Ref  17). 

When  present,  the  specifications  from  AFrDl,-TR -76-1  25  were  used 
since  this  report  specifically  addressed  RPVs.  The  handling  qualities 
in  the  Military  Speci ficat ions  were  used  to  complement  those  areas 
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which  ArFnL-TR-76-1 ?5  did  not  cover.  Table  IV  contains  the  results 
of  this  comparison  study.  It  consists  primarily  of  eigenvalue  specifica¬ 
tions  which  establish  the  allo'wable  range  of  values  for  the  roots  of 
the  characteristic  equation. 

Several  restrictions  in  determining  these  handling  qualities  were 
imposed  from  the  onset: 

--  The  general  specifications,  such  as  "Each  n.ode  will  have 
sufficient  damping",  are  not  repeated  in  Table  IV. 

--  Since  the  aircraft  in  this  study  was  a  RPV,  all  stick  force, 
pedal  force,  and  control  force  specifications  were  omitted,  it  was 
assumed  that  the  servo  v/ould  generate  sufficient  power  to  deflect  all 
control  surfaces. 

--  Since  this  RPV  was  designated  as  a  research  vehicle,  it  was 
assumed  that  it  would  be  flown  on  a  "good"  day.  Therefore,  atiDospheric 
disturbance  specifications  were  deleted. 

--  This  RPV  will  be  flown  by  sight;  thus,  display  specifications 
were  omitted. 

--  It  v/as  assumed  that  the  data  link  was  sufficiently  i'eliable 
and  capable  of  transsiitting  all  necessary  information. 

--  Redundancy  and  failure  rates  were  not  considered. 

Equations  of  Motion 

The  aircraft  equations  of  irotion  are  a  sot  of  nonlinear,  coupled 
force  and  moment  equations  which  completely  describe  the  motion  of  an 
aircraft  in  inertial  space.  By  assuming  perturbations  about  straight 
and  level  flight,  these  equations  can  be  decoupled  into  a  longitudinal 
sot  and  a  lateral -directional  set.  Linearizing  these  equations  about 
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Table  IV.  Hatidliny  Qualities 


Class!  fication:  Class  III  -  inediuin  [aaneuverabi  I  i  ty 


Flight  Phase:  Category  i\  -  rapid  ji’aneuvering 


Level:  Level  1  -  normal  operation 


Longitudinal  Case: 


''.sp 

t-nsp 

Residual  Oscillation 
Dynamic  Oscillation 


-  0.35  <  csp  <1.30 

-  See  Fig  5 

-  cp  >  0.04 

-  <  '  0.60  deg 

A  fhort  term  dynamic  oscillation 
produced  by  a  pulse  elevator 
shall  not  diverge  faster  than  a 
time  to  double  of  15  sec. 


!  Lateral -Di rectional  Case: 

'  ‘'DR 

1  '-^'nOR 

■  Dutch  Roll  Tiine  Parameter 

i  „ 

Spiral  Mode 

1 

) 

I  Roll  Effectiveness 

l 

Roll -Spiral  Coupling 


-  ’DR  >0.19 

-  'n[3R  >  0.40  rad/sec 

-  'DR-'-nDR  >  0.35  + 

where  =  0.014 

C-iiDR^  iT/fiDR  -  20] 

-  Tr  <  1.4  sec 

-  For  a  steady  state  disturbance 
in  of  20  deg,  the  time  to 
double  must  bo  greater  than 

20  sec. 

-  >  30  deg  ■>  in  1.5  sec 

-  None 


an  arbitrary  equilibrium  condition  porraitted  the  perturbation  ecjuations 
of  motion  to  bo  formed.  These  expiations  were  written  in  state  equation 
form  so  that  individual  trans.^'er  functions  and  modal  characteristics 
could  be  determined.  The  complete  derivation  of  the  perturbation  state 
equations  is  found  in  Appendix  3. 

M.jdal  Analysis 

The  models  of  Chapter  11  worai  used  in  the  perturbation  equations 
of  motion,  resulting  in  twelve  sets  of  state  equations,  six  for  each 
of  the  longitudinal  and  lateral -directional  cases.  Recall  Heat  a  :,t  ite 
equation  is  of  the  form 

X  =  A  X  +  8  u 

where  the  A  matrix  determ-’nos  the  characteristic  cqiwition  of  tlie  system, 
bach  A  matrix  was  used  in  a  rci'awi  por  (jrjgr.iim  (Ref  16)  which  de  teviii  ned 
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the  characteristic  equation,  and  hence,  the  r.odal  charactc'ristics  of’ 
each  system.  Table  V  summarizes  the  characteristic  equations  vdiile 
Table  VI  sunmarizes  the  modal  characteristics. 

A  coiiiparison  of  Table  VI  vn  th  Table  IV  indicated  that; 

--  'jp,  >'nDf^>  Tr  satisfied  their  respective  handling  qualities 
by  a  vnde  margin. 

--  Cp  satisfied  its  specification,  but  only  by  a  small  margin. 

--  cai^p  did  not  satisfy  its  handling  quality  in  cases  4  and  5,  and 
was  at  its  upper  li:,iit  in  the  other  cases. 

--  and  T3  did  not  satisfy  their  respective  spec i f icat ions . 

--  Since  all  the  longitudinal  roots  v/ere  stable,  no  dynamic 
oscillations  exceeded  a  time  to  double  of  15  sec. 

--  flo  roll -spiral  coupling  existed  since  each  mode  f’.ad  a  real 

root. 

--  Frequently,  residual  oscillations  are  tlv^  result  of  nonlinear¬ 
ities  in  the  control  system  (e.g.  dead  zone),  iience,  no  simple  way 
to  evaluate  these  effects  existed  within  the  scope  of  this  study. 

Two  handling  qualities  required  more  investigation  than  the  above 
comparison:  dutch  roll  time  parameter  ( tpR 'nrjfj)  and  roll  effectiveness. 
The  dutch  roll  time  parameter  handling  quality  is  given  in  terms  of  the 
phi  to  beta  ratio  ( j ..i/f  | qr)  .  This  ratio  may  be  found  by  determining 
the  magnitude  of  the  phi  and  beta  envelopes  at  an  instant  of  time  from 
the  rudder  pulse  time  responses  (Ref  4:673).  Utilizing  this  procedure 
and  the  time  responses  from  Appendix  C,  the  dutcii  roll  time  paraim.'ter 
specifications  were  determined  using  the  average  it/floR  with  the 
following  results: 
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Table  V.  Charac f;eristic  Ffjuations  of  the  Basic  RPV 


Longitudinal  Case; 

Case 

Characteristic  Equation 

1 

(s2  4-  6.34s  +  39.35)(s2  +  0.040s  t-  0.16)  =  0  I 

2 

(s2  +  7.38s  +  37.01)(s2  (■  o.042s  +  0.15)  =  0 

1 

3 

(s2  +  8.20s  +  24.16)(32  +  0.052s  +  0.11)  =  0 

4 

(s2  +  9.64s  +  65.72){s2  +  0.040s  +  0.18)  =  0 

5 

(s2  +  12.73s  +  31.49)(s2  +  0.046s  +  0.16)  =  0 

6 

(s  ^  11.27)(s  +  4.68)(s2  +  0.074s  +  0.093)  =  0 

1 

1  Lateral -Directional  Case; 

1 

C^se 

Characteristic  Equation 

1 

(s  +  11.50)(3  -  0.042)(s2  +  1.49s  +  19.03)  =  0 

2 

(s  +  11.18)(s  -  0.042)(s2  +  1.42s  +  17.76)  =  0 

3 

(s  +  10.88)(s  -  0.042)(s2  +  1.34s  +  16.50)  =  0 

4 

(s  +  11.30)(s  -  0.044){s2  +  1.23s  +  17.85)  =  0 

5 

(s  +  10.99)(s  -  0.043)(s2  +  1.17s  +  16.63)  =  0 

6 

(s  +  10.71)(s  -  0.044)(s^  +  1.08s  +  15.41)  =  0 
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Table  VI.  Modal  Characteristics  of  the  Basic  RPV 


1.63 
1 .87 


t 

t 

! 


t 


Case 

2 


Ihincn  iiK]  Qua!  i  ty 
'-DR >  0.48 
■"-DR'-nDR  >  0-78 


Comparing  these  values  \jith  Table  VI  shov/ed  that  the  dutch  roll  time 
parameter  specification  was  satisfied. 

Likewise,  tlie  roll  effectiveness  of  the  RPV  needed  to  be  determined. 
With  maximum  deflection  of  the  ailerons  (15  deg),  the  RPV  produced  62 
deg  and  64  deg  of  bank  in  1.5  sec  for  cases  2  and  5,  respectively  (Fig  6). 
Hence,  the  RPV  had  over  twice  the  required  roll  effectiveness. 

Although  not  specifically  mentioned  in  the  handling  qualities,  one 
other  characteristic  needed  attention,  the  ratio  of  the  oscillatory  com¬ 
ponent  to  the  average  ccnponent  of  bank  angle  ( losc^^ave^  •  ratio 

is  determined  by  analysing  the  peaks  of  a  bank  angle  time  response  due 
to  an  impulse  aileron  cor^nar.d.  Since  Cqr  <  0.20,  the  following  empirical 
formula  could  be  used: 


Jpsc  ^  ^'*'2 

'ave  .i;']  +'*^3 


where  fq  ,  ‘,2^  and  Pq  corresponded  to  the  first  second,  and  third  peaks 
(Ref  xxii).  Hence,  she  impulse  cime  responses  (Fig  7)  were  used  in 
Eq  (3)  and  produced 

Ca^e_  tosc/fave 

2  -0.0081 

5  -0.0108 


Thus,  the  RPV  was  not  very  oscillatory. 
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Ij'ansf^r  ^jnc_U_ons_  and  Tjnie  Responses 

ILefore  transfer  functions  and  time  responses  could  be  generated, 
the  actuator  dynamics  must  be  included  in  the  equations  of  motion. 
Modified  Futaba  FPS-14  servos  v/ere  used  for  the  RPV's  actuators.  Test 
of  these  servos  produced  two  different  representations: 


TP  -  (15^7)2 _  _  ^ 

servo  i2'+  (-21(781(1 5. 7)s  +Tr5.7*)2 

(14  45)2 

TFservo  -  ■s^~('2yr.8y(l4.45)s  +  (14. 4^) 2 


(4) 

(5) 


Equation  (4)  corresponded  to  a  5.5  in. -lb.  load  and  5.0  volt  input. 
Equation  (5)  corresponded  to  a  11.8  in. -lb.  load  and  10.0  volt  input. 
Since  the  RPV  would  have  a  peak  input  of  22  volts,  some  type  of  extra¬ 
polation  of  Eqs  (4)  and  (5)  was  needed.  It  was  decided  that  a  60  percent 
reduction  in  break  frequency  in  Eq  (5)  would  adequately  compensate  for 
the  larger  input  and  whatever  slippage  that  occurred  due  to  the  linkage. 
Hence,  the  following  representation  was  utilized  for  the  RPV  actuators: 


TF 


(7.22)2_  _ ^ 

servo  -  s7-+  l2Tr.WT.22ls  +TJ72W 


^2.20 _ _ 

sT  +  Tl.5'6  s  t-  52.20 


(6) 


The  correspondi ng  differential  equaticn  is 


5'=  -5--.?0i  -  ll.Sfj^  +  S2.20a,.n^d  (7) 

The  equations  of  motion  were  augiuontcd  with  Eq  (7)  resulting  in  the 
open  loop  state  equations  given  in  Fig  8.  From  tliese  equations,  the 
transfer  functions  for  thu  basic  RPV  were  geni?ratod  (Ref  16).  Only  Case 
2  is  presented  in  this  report  (T.'blc  VII)  since  the  FCS  design  will  I'O 
based  on  this  case. 
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-0.059  0.33  -0.29 


Table  VII.  Tran';for  functions  lor  the  basic  P.fV,  Case  2 
Lomji  tuciinal  Case: 

=  (s2  +  7.33  s  +  37.01)(s2  +  0.042  a  *  0.15)(s--  +  11. CG  s  ->  32.20) 

=  -2.34  (s  +  4.76)(s  -  163.08) 

=  -19.04  (s  +  114.53)(s2  +  0.57  s  +  0.28) 

=  -2219.32  (s  +  0.C89)(s  +  2.82) 

=  -2219.32  s(s  +  0.089)(s  4-  2.82) 

=  52.20  (s2  +  7.38  s  +  37.01)(s2  +  0.042  s  +  0.15) 

=  52.20  s(s2  f  7.33  s  +  37.01)(s2  +  0.042  s  +  0.15) 

Lateral -Directional  Case: 

V  =  (s  +  11.13)(s  -  0.042)(s2  +  1.42  s  +  17.76)(s2  +  11.56  s  +  -32.20)2 

N,®c.iid  " 

-  1.23  (s  -  2.27)(s  +  2.54)(s  -  26.33) 

^cmd 

^■acmd  "  t  1.43  s  +  15.19) 

^  '■  I  5  -I-  15.19) 

-  52.20  (s  +  11.18)(s  -  0.042)(s2  +  1.42  s  -i-  17.76)(s2  ^  II. -6  s  ^  ^^;'■;) 

'■^cmd 

il'Jc,i,d  "  ^2.20  3(s  +11.18)(s  -  0.042)(s2  +  1.42  s  +  17.7C)('.>’  *  II. 5C  v  .32. 2n) 

:r^»'  -  0 

'^^cmd 

n:'-  ^  0 

'acnid 

=  0.51  (s  -  0.078)(s  H  11.91)(s  79.23) 

■-  -40.04  (s  11.65)(s2  +  0.24  s  0.49) 

■  i^cnid 

N.-rcmd  "  (s  h-  4.32)(s  -  14.21) 

NPrcnid  =  +  4.32)(s  -  14.21) 

^*  rcmd  "  ° 

"Orcnd  ’» 

N.J^cind  ^  52.20  (s  +  11.13)(s  -  0.042)(s2  +  ^  /\2  s  +•  17.76)(s2  +  ll.'Si  t  '.."^.fO) 

■=  52.20  s(s  11.18)(s  -  0.042)(52  ^  1.42  s  17.7())(s''  ;  11. 5G  s  ^  '-.'.CO) 


n 


N 


f^cmd 

a 

‘I'ecmd 

9 

'?ecn,cl 
^  ^cmd 

f^e 

^  *^cmd 


o’" 

fcind 


Typically,  a  pilot's  coivicand  consists  of  an  input  and,  a  short 
time  later,  removal  of  that  input.  With  this  in  mind,  a  1  see  pulse 
was  selected  as  the  input  for  all  time  responses.  The  r:aijtiitude  of 
this  p(jlse  v/as  -5  deg  for  tfie  elevator  coii.mand  and  f5  deg  for  the 
aileron  and  rudder  coiiiinands.  A  set  of  tiii'e  responses,  for  case  were 
generated  (Ref  16)  using  these  inputs  and  are  contained  in  Appendix  C. 

Even  though  only  one  input  was  used,  several  general  observations 
were  made. 

--  All  the  time  responses  had  the  shape  that  would  be  expected 
:'ar  this  input  (Refs  2:47,  12'1-130  and  12:306,  668). 

--  The  peak  values  were  high.  This  stoim'.ed  from  the  slo\:  airspeed 
and  l->rge  control  surface  effectiveness. 

--  The  phugoid  nalu.-al  frequency  was  i:.uch  higlier  than  most 
di rcraf t . 

--  The  dutch  roll  n:o'le  was  very  oscillatory. 

--  The  spiral  "(jcie  had  a  dominant  effect  tho  t  could  i:ot  bo  ignored. 
C.)ric  ]  i.s  ions 

St7(':'al  cGnclusions  'wore  drawn  fio;:  t!ie  riodal  analysis  and  time 
resprrr.-'m  Tr 'Se  rneclusi.tns  define  t‘)e  u'osion  obje-Lti ves  for  Chapter 
IV.  In  ti’, ;  1'  i  ji  t'jfl'ino!  c'se,  r-nugni.l  r.alural  r>a;quoncy  needed  to  be 
dccrca  ,od  '’incii  o‘,c  i  1  !  a  t  ions  of  this  am.on  i  tud"  and  frequeiicy  were 
undesirable.  In  tiie  1  a  teral  -cl  i  rec  t  i  or.al  case,  tiio  effect  of  dutch 
roll  esr  illation  noed.d  to  be  reduced  and  the  spiral  i!:ode  needed  to 
be  stabi  I  i  .■’ed . 

Tile  effi'Ct  of  the  TVU  on  the  RPV  was  miuiii’al.  All  nal-ural 
freguL-nc,  i  ls  w"r(.-  slightly  reduced  as  was  the  short  period,  da'iiing 
ratio.  ihi-  f  ,er  ^diarac  ter  i  sties  remained  eisentially  unehinged. 
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A:,  xcn  location  inovtd  further  uft,  Lho  following  ciiarac ter i sti cs 
were  increasoh;  -sp,  •p,dnd,nQj^.  Also,  .n^pWas  decrea>ed.  All 
other  characteristics  rei::ai!iod  essc-ntially  unchanged. 
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IV.  fi^sinn  and  Eval nation 

An  dnsign  processes  consist  of  the  sa^:e  procedure:  synthesize, 
evaluate,  and  resynthesize.  The  purpose  of  this  chapter  is  to  complete 
tiiis  design  process  for  a  PCS  to  correct  the  deficiencies  found  in 
Chapter  HI.  Tlie  synthesis  phase  consisted  of  a  longitudinal  design 
and  a  lateral -directional  design.  Tiie  evaluation  phase  consisted  of 
three  sensitivity  studies,  a  final  handling  qualities  chock,  the  gen¬ 
eration  of  transfer  functions  and  time  responses,  and  a  nonlinear 
simulation. 

Long i tudi nal  Desj^n 

Even  though  the  longitudinal  cast-  satisfied  all  its  handling  qual¬ 
ities,  the  phugoid  :xde  v.as  quite  oscillatory.  Thus,  it  was  desirable 
to  reduce  this  esc i 1  la tory  response  eitner  by  increasing  p  or  de¬ 
creasing  ;.np.  Because  of  tne  constraint  of  only  rates  for  acceptable 
feedbacks,  the  choices  of  feodl)acks  for  aucii.eiitation  v.'cro  pi  ti.h  rate 
ana  picen  c,!;/,.  v/it'n  coi^pmsation. 

It  •  "jir-jd  .rily  a  i'ricf  glance  at  the  root  locus  of  tr.e  pitch 
rare  feC'Jbc.cr  cysr.'M  (Fig  9}  to  r-^alize  thal;  it  would  not  provide 
ap,)ropri  s  ce  cngi  .::,ldtii;n.  This  was  due  to  the  zero  that  \uis  very  close 
to  the  piingoid  unins  whic.i  tren:  I  a  tc.i  into  .a  7a*'y  large  gain  for  a  siiall 
inciease  in  campinc  ''atio.  In  itS'-lf,  the  large  gain  was  not  a  cwoblei::; 
rather,  it  was  its  efrext  on  the  other  roots  that  presented  th'.'  diffi¬ 
culty.  A  large  gain  v-anld  cause  the  short  period  root  to  becoi  r  Is'ss 
stable  and  its  hw'dling  quality  te  bGcome  unsatisfied.  In  addition, 
the  acturilor  dynumir.s  'would  chr.ngi>  greatlv  wliic‘"  W''-,  n-nies r  ib  1  e . 
Therefore',  tne  cnnclusic.n  v.as  t'nt  ths'  additie--  of  '•■ei : ',fi"m. ‘.’''u  would 
be  advan ta genus . 
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The  dosire  to  koop  tke  coi"pon';aticn  oii  plo  and  to  f?1i:;!in3to  the 
effect  of  the  zero  at  the  origin  dictated  the  ciioice  of  a  pure  integrator 
as  the  feedback  con:ponsator .  The  root  locus  (Fig  10)  revealed  that  the 
zero  at  the  origin  was  effer.ti vely  relieved  and  that  the  short  period 
root  b  !canie  unstable  faster  than  without  coi''i!)ensatiGn.  The  short 
period  handling  equalities  restricted  the  aihount  of  feedback  gain  to  less 
than  5.14.  A  feedback  gain  of  0.050  produced  the  best  cofspros.i so  between 
iiiaeiiiiL'tn  augv.ontation  and  minimuM  eff.-Ct  on  the  other  roots.  Thus,  the 
control  law  was 

"f^cnid  (s)  “  cin  (s)  f  ■  q  (s)  (B) 

An  exanination  of  ^-he  piuigoid  characteristics  showed  th.it  the  damping 
ratio  increased  frost  .-.p  -  0.055  to  op  =  0.25  while  the  natural  frequency 
remained  essentially  constant.  Hence,  the  oscillatory  phugoid  mode  was 
reduced  substantial ly. 

L  a  ttn’a  1  41  i  r  ec  tj  o  na  1  pe_s_i^i 

The  analysis  of  Chapter  III  suggested  that  the  dutch  roll  damping 
ratio  should  be  increased  and  the  spiral  mode  should  be  stabilized.  A 
s,;andjrd  yaw  damper,  cc.nsistiag  of  yaw  rate  feedback  to  the  rudder 
(passing  through  a  washout  circuit),  was  used  to  augment  the  dutch  roll 
mode.  The  main  difficulty  in  this  design  was  selecting  the  wasliout  pole 
location.  A  large  time  constant  for  this  pole  indicates  a  large  in¬ 
crease  in  both  damning  ratio  and  uncoordinating  effect  (Kef  12:f''T). 

A  washout  pole  location  of  s  =  -1.75  (T  ^  0.57  sec)  produced  tiie  l.r.st 
con;promise  to  the  trade-off  described  above.  A  feedback  gain  of  0.033 
was  selected  because  it  corresponded  to  I  he  r-aximum  dutch  roll  dansving 
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ratio  (see  Fig  11).  Thua,  tho  cor.trol  law  was 

I'Cind  ^in  ■*■  5  +  1  ’75 

An  examination  of  the  dutch  roll  characteristics  revealed  li.at  tlie 
damping  ratio  increased  fro.n  -  0.17  to  -  0.20  while  the  natural 
frequency  increased  from  ~  A. 21  rad/sec  to  ‘■n[)|^  ==  ‘-3. 12  rad/sec. 

To  siisplify  the  syste!ii,  tho  vrashout  circuit  was  rciiioved  from 
F  ,3).  1‘,'itn  this  control  lew,  a  focdhav'k  gain  of  0.032  produced  the 
r  .  u;n  dutch  roll  damping  ratio  (see  Fig  12).  The  dutch  roll  cliar- 

a... eristics  revealed  that  the  da,.;ping  ratio  increased  to  =  0.21. 
Untortunate-ly,  this  increase  was  not  enough  to  satisfy  the  handling 
qualities  by  an  adequate  margin.  Thus,  the  washout  circuit  was  rein¬ 
serted  into  Fq  (9). 

The  spiral  /rede  was  most  dominant  in  D/e  lateral  tine  responses. 
Thus,  either  integral  of  roll  rate  cr  yaw  rate  feedback  to  the  ailerons 
v/ould  be  the  liiost  effective  means  of  stabi  1  liring  tho  spiral  mode.  For 
the  integral  of  roll  rate  feeback  systeni  (Fig  13),  a  feedback  gain  of 
0.012  brcugfit  the  spiral  root  just  inside  the  left  hand  piano  resulting 
in  a  control  law 

‘•cmd  (^)  =  A>in  (s)  +  -----  P  (s)  (10) 

Implenientation  of  Fq  (10)  had  no  effect  on  any  1a teral -di rectiona  1 
characteristic  exc/'pt  the  spiral  root  wiiicfi  moved  to  s  =  -0.0003  fro.:) 
s  -  0.041 . 

For  the  yaw  rate  feedback  system  (Fig  14),  a  feedback  gain  of 
0.044  ostablisiiod  a  stable  spiral  i.,ode. 
system  was 
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Tho  control  law  for  this 


Figure  11.  Root  Locus  for  the  Yaw  Damper  Including  the  Washout  Circuit 


(.1 


Figure  12.  Root  Locus  for  the  Yaw  Damper  excluding  the  W,:  diout  Circuit 
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Figure  1^-.  Root  Locus  for  Yav;  Rate  pLcdback 
to  the  Ailerons 
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^^cmd  (s)  -  ^in  (s)  -  0.044  r  (s) 


(11) 


Implementation  of  Eq  (11)  also  had  no  effect  on  any  1 atcral -di rectional 
characteristic  other  than  moving  the  spiral  root  to  s  =  -0.0001  from 
s  =  0.041. 

Thus,  both  Eqs  (10)  and  (11)  accomplished  the  desired  goal  of 
stabilizing  the  spiral  mode.  Equation  (11)  v/as  selected  as  the  more 
desirable  control  law  because  it  did  not  require  any  compensation  or 
any  additional  instrumentation,  such  as  a  roll  rate  gyro. 

Sensi ti vi ty  Studies  and  Redesign 

Three  types  of  sensitivity  studies  were  used  to  evaluate  the 
PCS:  a  flight  condition  study,  a  parameter  study,  and  a  physical 
quantity  study.  Each  study  based  its  results  on  a  modal  analysis. 

The  flight  condition  study  varied  the  airspeed  throughout  the  entire 
speed  regime.  The  parameter  study  varied  each  parameter  in  the 
equations  of  motion  (except  Uo>  q»  S,  and  c)  by  25  percent.  Lastly, 
the  physical  quantity  study  varied  the  vertical  tail  area. 

The  modal  analysis  from  the  flight  condition  study  showed  that 
three  handling  qualities  were  not  satisfied  over  the  entire  airspeed 
regime.  First,  the  short  period  natural  frequency  handling  quality 
became  unsatisfied  at  the  higher  airspeeds  (M  i  0.12).  Feedback  gain 
adjustment  had  no  effect  on  increasing  the  range  of  airspeeds  that 
satisfied  this  handling  quality.  Titus,  a  restriction  was  iinpos'^d 
stipulating  that  this  handling  quality  vtould  not  be  satisfied  at  the 
higher  airspeeds.  This  restriction  was  acceptable  when  considering 
that  the  higher  airspeeds  were  clearly  in  the  upper  limits  of  the 
RPV's  cartabi  1  i  ties .  Secondly,  this  study  revealed  that  the  spiral 
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mode  wa ;  too  unstablo  at  low  airs|)Cods.  InortMO i ny  th'.‘  foodback  gain 
in  Eq  (11)  to  0.09  stabilized  this  mode  throagi-oj t  the  cntifo  airspeed 
regime.  Lastly,  this  study  found  that  the  dutch  roll  damping  ratio 
did  not  satisfy  its  specification  at  the  higher  airspeeds  (M  b  0.15). 

A  feedback  gain  of  0.03  maximized  this  dampinn  ratio  throughout  the  air¬ 
speed  regime.  Unfortunately,  as  in  sfiort  f'Oriod  natural  frequency,  the 
entire  airspeed  rogiire  could  not  be  satisfied  and  tlie  higfter  airspeeds 
restriction  was  ia. posed.  Again,  this  was  acceptable  when  consic.'ering 
the  RPV ' s  capab il i ties . 

Hence,  the  control  laws  (E|S  (3),  (9),  and  (11))  were  refined  to 

^eemd  (s)  =  Ejn  (s)  +  <1  (s)  (12) 

l^cind  ('^i  ‘  i^in  (s)  -  0.090  r  (s)  (13) 

'-if'ciiid  (s)  =  Rin  (5}  ^  (1'^) 

Plots  of  the  pole  locations  due  to  a  varying  airspeed,  with  the  above 
control  laws  iniplemcnted,  may  be  found  in  Fig  15.  These  plots  show 
pole  locations  for  airspeeds  from  Mach  0.05  to  Mach  0.15  at  Mach  0.01 
intervals.  The  arrows  indicate  increasing  airspeed. 

The  second  study  entailed  tabulating  the  modes  of  the  RPV  as  each 
parameter  in  the  equations  of  motion  (except  Ug,  q,  S,  and  c)  was 
varied  .':25  percent.  Comparing  the  extremum  (Table  VIII)  with  the 
handling  qualities  (Table  IV)  showed  that  all  handling  qualities  were 
satisfied.  Thus,  an  error  in  any  one  parameter  of  less  than  25  percent 
will  not  cause  any  handling  quality  to  ho  exceeded.  Of  course,  no 
assurance  can  be  made  if  mote  than  one  parameter  changes. 
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Table  VIII.  Extremun  of  the  Parameter  Sciv-,i  ti vi  ty  Study 


Source  of  the  extremum  is  given  in  parentheses.  The  sign 
ifidi cates  whether  the  parameter  was  increased  or  decreased 


The  last  study  was  completed  because  of  the  single  vertical  tail 
approximation.  Three  scale  factors  (1.50,  1.75,  and  2.00)  were 
examined  to  determine  the  effect  of  the  vertical  tail  area  on  the  system 
modes.  As  expected,  the  vertical  tail  area  had  no  effect  on  the  long¬ 
itudinal  case.  Little  effect  was  found  on  tlie  lateral -directional  case. 
A  plot  of  pole  locations  due  to  varying  vertical  tail  area  is  shown  in 
Fig  16.  The  arrows  indicate  increasing  vertical  tail  area.  Comparing 
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Figure  16.  Pole  Movement  due  to  Variations  in  Vertical  Tail  Area 


Fig  16  and  Table  IV  revealod  that  all  handling  qualities  were  satisfied. 
Thus,  a  substantial  variotion  in  vertical  tail  area  diu'  rmt  cause  any 
handling  quality  to  becona  unsa tisf ivcd . 

ilandjing  Oual  i  ties  ChecJ^ 

A  final  check  of  tne  handling  qualities  was  ii.ade  with  the 
control  laws  (Eqs  (12)  -  (14)).  A  comparison  of  tfie  augmented  RPV, 
the  basic  RPV,  and  the  handlinij  qualities  (Tal)le  IX)  reveal(.“d  Lhai  all 
characteristics  '.atisfied  their  r.espcctive  handling  qualities. 
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Taljlo  IX.  Coii'parison  of  tho  Ai,'<:;;;cnf:od  RPV, 
Basic  RPV,  and  tho  Handling  Qualities 


Augmented 

Basic 

llandl  ing 

RPV 

RPV 

Qudl i ti os 

?sp 

0.58 

0.61 

0.35  c  fsp  <1.30 

'^’njp,  rad/sec 

5.74 

6.08 

Fig  5 

Cp 

0.25 

0.055 

Cp  >  0.04 

wnp,  rad/sec 

0.41 

0.39 

★ 

Tr,  sec 

0.087 

0.039 

Tr  <  1.40 

Ts,  sec 

22.91 

-24.04 

Ts  >  0  or  Ts  <  -23.85 

(^•) 

(16.66) 

(>20.00) 

^DR 

0.31 

0.17 

‘■■DR  >0-19 

(^•'noR,  rad/sec 

4.60 

4.21 

’''nQR  >  0.40 

Dutch  Roll  Time 
Parameter,  sec"^ 

1 .44 

0.71 

CDR  ‘^noR  >  0.35 

*  No  iiandling  quality  exists. 

**  Negative  sign  indicates  instability.  Time  to  double  (sec)  is  given 
in  parentheses. 
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As  before,  tlie  dut'ch  roll  ti!':e  p.irr:' eter  (  l,.;;  riijj|;/  '  j  ec  i  fi  sab  ion  was 
dei'endent  on  the  phi  to  beta  racio  (;;/;-1dr)-  ti'-s  responses  of 

Appendix  C  revealed  that  the  yrap'iical  method  of  d-;ter"ii ni 
in  Chapter  III  could  not  be  utilized.  Therefore,  an  analytical  i-ethod  vjas 
used.  The  |[jR  is  the  ratio  of  the  magnitude  of  the  ,  ccirponent  to 
the  magnitude  of  the  ;  component  of  the  dutcfi  roll  eige.ivec tor .  Thus, 


Since 


1  O.OuCT  +  j"O.Ob67 
=  0.13 


2.84 


I  20.00 


the  dutch  roll  time  paraiaoter  specification  became  'pR  '^0.33, 

which  Wci.s  satisfied  by  the  augmented  RPV. 

In  addition,  the  roll  effectivimioss  of  the  augmented  RPV  needed 
to  be  computed.  With  l•^'xi;llUl:l  deflection  of  the  ailerons  (15  deg),  the 
aufpi'.entcd  RPV  proilucsJ  62  deg  of  bank  in  1.5  soc  (Fig  17).  Hence,  no 
change  in  roll  uf fccti vtoness  was  experienced  due  to  the  FCS. 

Lastly,  since  gj.;  0.20,  Lq  (3)  could  not  be  utiliz  d  to  detotiiiine 
•‘osc/''ave'  Rather 


‘osc  _  *1  - 

'  .1  ve  ■  1  f  ■  2 


siioiild  be  used 
was  used  in  Lq 


(Ref  4:  xxii),  Ttur, ,  b.he  impulse 


(15)  and  p.rodijred 


•  <v,(: 

•  ave 


0.044 


ti.'-.e  rfs[;iwr.e 


(I’d  IB) 
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Kd'icp,  ti  0  au'jroiited  RPV  was  not  very  oscillatory. 


^  ‘'11'^  Time  Response's 

Before  transfer  functions  and  time  responses  could  be  determined, 
the  closed-loop  state  equation  must  be  derived.  Recall  the  form  of 
the  open-loop  state  equation 

“  A  x_  +  S  (16) 

where 

u  =  K  X  +  IN  (17) 

Substituting  £q  (17)  into  Lo  (16)  and  simplifying 

X  =  Ax  +  B(Kx  +  ]fl) 

=  (A  +  BK)x  +  B  IN  (18) 


Substituting  the  control  lows  (Cqs  (12)  -  (I'l)  into  iln  (18)  provided 
tiiC  closed-loop  state  equation  given  in  Fig  lt>.  Notice  that  a  new  state 
variable  (lad  to  ha  added  to  the  lateral -directional  state  eqiiation  due 
to  the  v/ashout  circuit  in  the  yav.'  damper.  The  new  state  was  cliosen 
to  be  t yaw  damper  foed.b.ack  C  and  defined  as 


!:(.  ]  -  0.030  s  ^ 

s  V  1.75  '  ^ 


Using  the  closed-loop  state  equation  (Fig  19),  transfer  functions 
and  time  responses  of  tne  aug;,''nt.od  RPV  va.'re  uc  nc-a  t' si  (Ref  lb).  Ta,t)lo 
X  contains  the  resultant  tran  .fi^r  functions  wtn'le  tbe  time  rrsponsos 
are  in  Appendix  C.  Co;a[?arir.g  t‘''e‘^e  tiise  respr'uses  with  Lht?  liasic  RPV 
time  re  .ponses  verified  t' at  the  li.-sign  I'l.i  j  .'c  ti  v>'s  p.aci  bn  ti  actiieved. 


Tabl(?  A.  Tra.’i f'j'f'  J’ufictioris  for  tl'C  Ai;- .:'on  tod  RPV 
Lonqi  ti:  J  i  ndl  Case; 


1  v 

0.21  5 

)  +  0. 1  7) f s-  -i-  6 . G2  s 

+  32.33)  (s'^  +  12 

=  -2.34 

(s  + 

l./6)(s  -  163.08) 

^'Ein 

=  -19.04 

^  (s  + 

114.58)(s2  +  0.57  s 

+  0.28) 

1  0 

; 

-  -2219. 

32  (s 

+  0.0S9)(s  +  2.32) 

'  :i9. 

tin 

=  -2219. 

32  s(s 

;  +  O.G89)(s  +  2.82) 

iCe 

=  52.20 

(s2  + 

7.38  s  +  37.01 )(s2  + 

0.G42  s  +  0.15) 

!  .■'fl 

1  Nse 

-  52.20 

s(s2  1 

-  7.33  s  +  37.01)(s2 

+  0.042  s  +  0.15) 

i 

Lateral -Di roctional  Case: 


-  (s  +  11.50)(s  +  0.04'Ofs  +  2.32)(s2  +  2.89  s  +  21.12) 

(s2  +  11.56  s  +  52.20).s2  +  9.11  s  +  31.49) 

-  -64.78(s  +  0.21 )(s  +  2.42)(s  -  36.72)(s2  +  10.95  s  +  46.42) 

=  64.09(s  +  l.75)(s  -  2.27)(s  t-  2.54)(s  -  26.4^1 
(s2  -I-  11.55  3  +  52.18) 

=  2283.02(5  +  2.3G)(s^  +  2.96  s  +  20.39)(s2  v  9.49  r>  +  29.63) 

-  2288. C2  s(s  +  2.30)(s2  +  2.95  s  +  20.39](s^’  +  9.49  s  +  29.63) 

-  52.20(s  -  O.G41)(s  +  2.31)(s  +  11.1G)(32  -  9.32  s  +  30.75) 

(s2  +  2,63  s  +  23.44) 

=  52.20  s(r,  -  0.041  )(s  +  2.31  )ls  +  11.10)(s2  +  9.82  s  +  30.75) 
(s^  +  2.68  s  +  23.44) 

=  100.27  (s  -  O.GOG4)(s  -  Z.27){s  +  2.55)(s  -  26.32) 


Table  X.  (Continaod) 


i  Ain  =-■  100.27  s(3  -  C.000-])(s  -  2.27)(3  +  2.l.T)(s  -  25.32) 

:  nL  -  1.'^-’  (s  -  Q.^:)O0)(s  -  2.27)(5  +  2.55)(3  -  2G.32) 
i  ’  (sO  r  11. o6  s  H  52.20) 

I  -  26.05  +  0.019)(3  -r  l./5)(s  +  n.74)(3  +  79.24) 

I  '  (s^  +  11.63  s  +  51  .16) 

Hr.  -  -2000.00  (s  +  1.75)\S  n.65)(32  +  o.24  s  +  0.49) 

(s-  +  11.56  s  +  52.20) 

N,^.^  =  713.53  (s  I-  1.75)(s  10.47)(s  -  12.95)(s2  +  4.15  s  +  22.30) 

nS.  =  713.53  3(s  +  1.75)(s  +  10.47)(s  -  12.05)(s2  +  4.15  3  +  22.30) 
'^1  n 

=  9313. S4  (s  +  1.75', (5  11.63){s2  +  0.24  3  +  0.49) 

j  n 

i  Hiia  =  9813.84  s(s  +  1.75)(s  +  11.55)(s2  +  0.24  s  ^  0.49) 

I  ^i'l 

1  Np*;'  52.23  (s  ^  0,045)(5  +  1.75)(s  +  11.52)(sl^  +  1.37  s  +  16.24) 

i  111  (s'-  +  11.19  s  -I-  53.26) 

ri'i”  -  52.2  s(s  H  0.04.5)(s  +  l.75)(s  +  ?1.52)(,-;2  +  1.37  s  +  16.24) 
'Mn  (s2  +  11.19  s  53.26) 

Nn.  -  -62.70  (s  -  0.0003)(s  +  11.55)(s2  +  0.24  s  i-  0.49) 

I  (sO  +  11.56  s  +  52.20) 
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Non!  ir;ear  '^ii’iuLilion 

The  last  fona  of  evciluati  jo  ■.•.as  a  tG;  ,pli_'te  nonl  inear  s  Ir/jlation 
which  examined  the  effect  of  cuai'Ming  on  the  fCS.  Initially,  an 
analog  simulation  was  cGinpletod  'with  the  r'‘sults  being  of  low  quality. 
Therefore,  a  digital  siiculation  was  accomplished  using  the  numerical 
integration  subroutine  ODE  (Ref  13)  (see  Fig  30).  The  equations  of 
motion,  in  their  raost  basic  form  (Eqs  (35)  and  (33)),  were  used  as  was 
the  same  control  inputs.  Recall  the  nonlinear  states  (except  U)  have 
the  same  value  as  the  perturbation  states  since  the  equilibt'ium  con¬ 
ditions  i’ere  chosen  zero.  Appendix  C  contains  the  results  of  this 
nonlinear  simulation.  A  comparison  with  the  other  tii',8  responses 
indicated  that  the  coupling  effect  between  modes  was  quite  small. 

Thus,  the  decoupling  assumption  produced  mininial  error  and  was  accept¬ 
able. 

In  the  foriiiulatlon  of  the  simulation  algoritlv.n  (Fig  20),  one 
critical  selection  was  made,  the  increment  time.  Shannon's  sampling 
theorem  dictated  that  at  d  0.68  sec  for  complete  reproduction  of  the 
response  (Ref  9:  Sec  III,  21).  In  practice,  a  much  faster  increment 
time  thai\  At  -  0.68  sec  is  desirable.  Thus,  At  -  .02  sec  was  chosen 
for  the  first  2.0  sec  of  the  simulation  and  At  -  0.10  sec  for  the 


remainder  cf  tlie  simulation. 
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V.  A  Practij:aj  -ApPj 

This  study  would  not  be  complete  without  adfiressing  an  application 
of  this  RPV.  It  was  hoped  that  by  demonstrating  its  usefulness,  fnotiva- 
tion  for  continued  research  in  this  area  would  bo  generated.  The  chosen 
application  was  torned  model  matching,  i.e.,  making  the  RPV  have  similar 
dynamic  characteristics  as  another  aircraft  through  the  design  of  an 
expanded  PCS.  The  aircraft  chosen  for  duplication  was  the  F-15.  The 
first  step  of  this  procedure  was  to  determine  the  F-I5's  dynamic  char¬ 
acteristics.  Then,  three  design  procedures  were  investigated  for 
accomplishing  this  duplication:  Guillemin-Truxal  design.  Entire 
Eigenstructuro  Assignment,  and  an  extension  of  Observer  theory.  Only 
the  longitudinal  case  was  examined  since  the  lateral -directional  cose 
would  be  a  repetition  of  the  same  theory. 

AnjJj/sij  0^  the  F-15 

The  mathematical  model  for  the  F-15  was  obtained  froiii  Ref  19  and  is 
contained  in  Table  XI.  As  with  the  RPV,  the  equations  of  motion  were 
augmented  with  the  actuator  dynamics.  The  following  "oprosentation 
was  utilized  for  the  F-15  servo  (Rof  19): 

—  20 
' '"servo  =■- 

The  corresoonding  differential  equation  was 

66  -20  '.e  '■  20  ''OQiiid  (19) 


Iho  equations  of  motion  were  augiiien tc-d  with  Fq  (19)  resulting  in  the 
longitudinal  state  equation  given  in  Fig  21.  Trarsfor  functions 


lablc!  Xl.  F-15  Matr.c-ii.dtical  Model 


ni  =  10B5  sluqs 
Uq  =  655  ft/sGC 
q-  =  338  lb/ft2 
S  =  608  ft2 
c  =  16  ft 

ly  -  155746  slu(j-ft2 
Cxu  =  --04 
Cv-  =  0 
Cx  "  0.24 

a 

Cxq  =  0 

Cx5e  “ 

=  -0.30 
Cl.  =  1.01 
Cl'‘  -  3.78 
CLq  ^3.12 
C2,5e  =  0 
C(.iu  =  0 
Cm-  =  -1.20 
Cmq  =  -3.70 
Cm.5e  "  -0-65 


stability  ond  control  derivatives  have  the  units  of 
rad . 
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(Table  XII)  and  tiiiio  responses  (Figs  25  -  29)  were  generated  frciii  this 
state  equation  (Ref  16).  Also,  the  longitudinal  characteristics  were 
determined  as 


i-risp  =  2.90  rad/sec 

?p  ==  0.10 

(■iOp  =  0.063  rad/sec 
Gui  1 1  enii n-Truxal  Design 

The  Guil  leiT.in-Truxal  design  procedure  is  basically  a  pole/zero 
placement  technique.  From  a  desired  closed-loop  transfer  function, 
the  appropriate  cascade  or  feedback  compensator  is  mathematically 
calculated  (Ref  6:408).  This  procedure  guarantees  duplication  of  the 
desired  closed-loop  transfer  function  for  single  input-cingla  output 
systems.  However,  in  this  application,  the  RPV  was  a  single  input 
(4ec„id)-niul  tiple  output  ('u,  a,  o,  q,  -:e)  system;  therefore,  the 
validity  of  this  technique  was  questionable.  However,  it  was  decided 
that  the  duplication  of  ono  state  iiiiglit  prove  valuable.  Of  course, 
the  duplication  of  one  state  v^ould  alter  the  other  three  states  which 
might  provide  a  closer  match  to  all  of  the  F-15  states. 

For  a  cascade  compensator  Gc(s),  the  closed-loop  unity  feedback 
transfer  function  is 

M(s)  _  Gr(s)  G(5) 
i)(s)  “  1  +  (3c(s)'G(sT 

Solving  for  Gc(s) 

r  ( <•)  -  ( 5 ) 

■  [D(s)  -  N(s)1G(sT 
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Applying  this  equation  to  the  RPV  and  F-lb  state,  (Rof  16)  r-;sultc(l  in 
four  cascade  compensators  (Table  XII!),  all  of  which  wen;  very  compli¬ 
cated.  From  cs  hardware  point  of  view,  these  compensa tors  would  bo 
extremely  difficult  to  impl e;r:cnt ,  if  not  impossible.  Thus,  cascade 
compensators  were  rejected. 

This  sa:".e  prscedure  was  repeated  for  feedback  coi.  ponsa  tors .  The 
closed-loop  transfer  function  is 


ri(s)  _  Gjs) 

b(s)  ■  1  + 


Solving  for  Gc(s) 


F‘C  ( s ) 


l)(s)G(s)  -_N(s_l 

n(syG'(s}' 


This  equation  resulted  in  four  feedback  cempensators  (Table  XtV)  wliich 
were  simpler  than  the  cascade  compensators,  but  still  impractic;al  to 
implement.  Since  botli  coicpensa tors  were  not  feasible,  too  Guillemin- 
Truxal  design  prncodtire  was  dismissed  as  indpp»-opriato  for  this  3i)[)li- 
cation. 


En_Uf2e  Fi  gen  structure  Assignrnonjt 

Fntire  Fi  gens  true  ture  Assignment  is  a  des  1  auc  wherc'by 

the  engineer  may  assign  both  the  eigenvalues  and  c ig-'^nvf'c ton,  to 
system.  This  technique  is  bas^d  on  full  state  fetr'f'ark  an  )  u'  c;  ,  the 
dtsired  eigenvalues  and  eig'’nv(;ctor ,  to  calculate  the  state  f' •  t  t 
matrix  (Pef  10).  With  cnl/  one  input  as  in  this  <!p;'l  ica  t  i  on ,  t'  ■ 
selection  of  the  closed- i.iop  oigc-nvalufjs  coisplePoly  deter'  ino.  tia’ 
eigi^nvectors.  Because  of  tin's,  the  effectiveness  of  tins  te;hnii]:r- 
was  also  questionable. 
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Pt'‘' v-fi  tl y,  t.!,'?  RPV  .ifid  the  F-15  are  different  ordc-r  systenis  duo  to 
tre  aetii.ilnr  reiu'es.n'.tationu  of  each.  Thus,  the  use  or  this  technique 
will  he  'iitfie  lit.  To  sii'.plify  the  process,  it  was  decided  to  employ 
a  reduced-oni-r  >  n.ud.'M  for  thr?  RPV.  The  criteria  used  for  the 

selection  of  tm.'  r 'sducoa  order  model  v;os  time  response  to  a  pulse  input 
The  Output  of  tne  three  , models  listed  below  was  plotted  against  the 
KFV's  second^ordsr  model  (Fig  ?2)  (Eg  (6)). 


Mo_de_l_ 

1 

2 

3 


R educed-0 rder  R c pres entation 

7.53_  _ 
s'  +  7.53 


S  +  6 
5 

S'+  5 


6  •) 


CJ 


/f?  ! 

'i!  i 

:  ,  / 


Model  2 


Model  3 


/ 

0 


Model  1  ' 


Second-Order  Model 


0.4 


1  .6 


0.8  1.2 

Time  (sec) 

Figure  22.  Reducvd  Order  Servo  A.ppro.-.imation  for  the  R!’V 
lodel  1  lias  the  sai.e  2  pcrMor, t  'U-ttine  tiiso  as  the  second-ordr-r  model. 


TiiO  other  eodnl,,  botli  of  vliich  h..id  slowm'  responses,  were  examined 


because  they  ;;3uld  be  isorc  likely  tu  pi'oduc^:  better  approxi.  .rition  to  th 
scccnd-ordor  model.  It  v/as  decided  tlict  moiel  2  provided  tl:e  icost  ac¬ 
curate  approxiriation  to  the  second-order  ;.:oo'cm  ,  so  it  was  used  thrcurjhout 
the  rest  of  tliis  application. 

The  F-15  eigenvalues  v/ere  assigned  as  the  closed-loop  eigenvalues 
of  the  RPV  (Ref  3).  As  expected,  the  associated  eigcnvectof’s  did  not 
correspond  to  the  F-15  eigenvectors  (see  Table  XV).  Thus  the  L'ntire 
r.i gens tructure  Assignment  design  technique  was  also  disisissed  as  inap¬ 
propriate  for  tin's  application. 

[^tension  cf_  Observer  Theory 

Witii  both  the  Gui  1 1  einin-Truxal  and  the  Entire  Eigenstructure  Assign¬ 
ment  techniques  yielding  unacceptable  results,  scsiething  more  basic  to 
control  system  design  was  needed.  Observer  t'neory  gave  the  necessary 
insiglit.  The  generalised  block  diagram  (Fig  23)  for  this  techriique  is 
in  terms  of  an  actual  model,  a  desired  piauel ,  and  several  reedback 
matrices  (L,  K,  and  K').  In  an  observer,  K  and  K'  would  be  identity 
matrices  and  tlie  error  feedback  v.ould  go  to  the  ciesirod  r:odel  instead 
of  the  actual  model,  i.e.,  the  desired  model  \.'ould  ho  rndo  to  track 
the  actual  nodel  .  in  isodel  iiK.tching  the  goal  is  for  tb,  '  actral  model 
to  crack  the  desired  piodel ;  he.nco,  the  error  feocioack  geos  to  the  actual 
m.odel  .  K  and  K'  wore  added  to  she  system  to  increase  flexi  Ir' 1  i  tv .  In 
this  apj.' 1  ica ti on ,  the  RPV  woulfi  be  the  actual  mO(iol  and  tr:e  P  -l  t  mrulcl 
be  the  desired  model. 

Thi:  clo'ed-loop  syst.^m  v.'as  expressed  in  forms  of  several  cem,  uf 
statii  variables.  It  v,  i  s  found  tint  the  best  choice  of  stale',  v/ru'i:  the 
five  desirr-d  aodel  (F-IS)  state  variables  (x')  and  th:}  five  errer  state 
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Figure  23.  Generalized  Bloc!;  Diagram  for  Model  Matching 
variables  (c),  v/here 

£  =  X  -  x'  (20) 

Thus,  the  closed-loop  state  equation  is 


rA  -  BLKC  I  (A  -  BLKC)  -  (A'  -  BLK'C) 


'bGc  +  B'" 

+  -  u  (21) 

B' 


The  servo  representations  have  been  included  in  A  and  A',  tjso, 
for  single  input  systems,  B  and  B'  will  always  bo  zero  vectors  except 
for  the  last  elcinont.  Hence,  K'  added  no  additional  floxibilily  to  the 
system  and  v.'as  assun^ed  to  l)e  the  identity  matrix.  For  this  application 


C  and  C  v;ere  also  identity  i;;a trices.  Final]/,  defining  R  -  LK  pro¬ 
vided  insight  into  choosing  the  fe-dback  matrices.  Tfius,  Fq  (21) 
became 


A  -  BR  1  (A  -  BR)  -  (A'  -  BL) 

1  _ 

r 

_ 

X  ' 

L-J 

r 

0  1  A' 

+ 


(22) 


This  form  of  the  state  equation  yields  insight  into  the  proper 
choice  of  weighting  matrices  (R,  L,  and  Gc).  The  solution  of  Eq  (26) 
is  in  terms  of  the  eigenfunctions  (vj  Thus,  R  should  be  chosen 

such  that  the  domir,  .nt  terms  of  tlie  eigenfunctions  of  (A  -  BR)  dio  out 
rapidly.  Eviari  though  all  eiger.valuas  of  (A  -  BR)  need  to  be  negative, 
this  does  not  indicate  tnat  all  eigenvalues  need  to  be  far  in  the  left 
hand  ijiane.  Only  the  ones  associated  with  the  dominant  eigenfunctions 
need  liave  large  negati/e  values.  [(A  -  BR)  -  (A'  -  BL)]  is  a  coupling 
matrix.  Ideally,  the  best  choice  for  this  matrix  would  be  a  zero 
matrix.  Unfortunately,  L  dees  not  have  enough  va'^iablcs  to  accomplish 
this.  Therefore,  L  should  be  ciiosen  to  iiiinii:ii.ze  the  coupling  effect. 
Finally,  it  would  bi:  desirable  if  an  input  did  not  directly  affect 
the  error.  Thus,  should  bo  crosen  to  produce  a  zero  vector  for 
(BGc  -  B'). 

The  techniques  for  choosing  R  and  G^  are  relatively  simple. 
However,  minimization  of  the  coupling  effect  represents  a  greater 
challenge.  To  simplify  the  notation  of  this  derivation,  define 
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=  A  -  BR 


b  =  (A  -  BR)  -  (A'  -  BL) 


Thus,  Eq  (?.2)  becomes 


£ 

a  ,  B 

_ 1 _ 

1 

V  ' 

L- . 

0  [  A' 

L-  J 


i  BGc  -  B' 


B' 


(23) 


Examining  the  zero  input  response  to  Eq  (23)  produced  valuable  insight 
into  the  importance  of  f.  The  zero  input  response  is 


i^''4) 


r(t) 

10 

'it  j 

"  l(0)' 

= 

I 

e  yj  fi 

x'(t) 

i  =  l 

x'(0) 

where 


>.i  ■-=  i  th  eigenvalue  of 
v-j  =  i  th  eigenvector  of 

r-j  =  i  reciprocal  basis  vector  of 


Ct  1 

B 

1 

— 1 
-  ?  1 

A' 

'  1 
_ l_ 

B  ‘ 

r 

0  1 
. —  1 

A' 

J 

Because  of  the  upper  block  Form  of 


a  ,  b 
0  1  A' 


0 


;  A- 

I  ^ 


]  ,  the  syst(:!n  eigenvalues  are 


simply  the  eigenvalues  of  ■<  and  A'.  Also,  Vj  r-j"*”  is  a  10  x  10  matrix 
v/iiosG  upper  right  5x5  block  (defined  as  ■ )  needs  to  be  minir.ized  for 
all  i. 

For  i.he  eigenvalues  coi^respondi rig  to  't,  the  eigenvreto'^s  and 
reciprocal  basis  vc-ctors  are 
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(25) 


n 

ii 

Li  “ 

0 

wi^D[Ai^  I  -  A-]-1 

where 

=  ith  eigenvector  of  a 
wj  =  reciprocal  basis  vector  of  a 


For  the  eigenvalues  corresponding  to  A',  the  eigenvectors  and 
reciprocal  basis  vectors  are 


y-i  = 

['MA'  ^  ^Aii' 

“I 

I 

ri  = 

0 

J 

W-j  ' 

(26) 


Where 

ij'  =  i^*^  eigenvector  of  A* 

wi '  =  i^^  reciprocal  basis  vector  of  A' 


Now  that  the  eigenvectors  and  reciprocal  basis  vectors  have  been 
determined,  y  needs  to  be  examined.  Recall  that  L  yields  control  over 
only  the  fiftii  row  of  .  Therefore,  the  equations  were  simplified  to 
isolate  these  terms.  For  the  eigenvalues  corresponding  to  u 


r  - . 


Vj  r 


i^.i 


.T 


[wiT  ;  .  A-j-l] 


or 


"  T,i  wi^  A  j  -  A-]-'  ^  (P/) 

Defining  -  [  '  i  ^  I  -  A']'^  and  expanding  Fq  (?7)  for  the  r-c  component 
of  y 


/O 


^rc 


k^l  '•If' 


i  r  'f'i  k 


5 


'kj  'I'jcj 


■■  j^=i  ^ir  wik 
or,  in  matrix  form 


5 

'i-jc 


■*■  f  i  r'^'f  i  3 


^5j  Pjc: 


Yrc  “  [firWiShc  •  •  •  fir^'fiSfScJ 


.^T'if'^ik  jj  =  i 


For  the  eigenvalues  corresponding  to  A' 


■51 


'■55 


Vili 


.T  _ 


[Ai/\ '  I  -  u]  ^  SCi ' 


i-.i 


Co"^!  Wi'Tj 


(28) 


or 


Y  =  [XiA'  I  -  u]'^n  ii’wi'T  (29) 

Defining  =  [XiA'  I  -  expanding  Eq  (29)  for  the  r-c  component 

of  Y 


5  5 


^rc  “ 

kM  j  =  l 

i'rj«jl<]  Oik’ 

wic ' 

r-. 

5  4 

k^l  ^  j^l 

'XrjjjkJ  "ik' 

5 

''■ic'  v'r5>  5k 

■■  i  k  '  i  c 

or,  in  matrix 

form 

‘''51 

r  •  r5 

•  re  L  ' 

ij'wic  •  •  • 

’  r5'  i 5 ' -''ic] 

1  ‘  55 

5  4 

1  f-jX]  'rj.  jk]  '.ik'^ic'  (80) 


/I 


Augnontiriq  Eq  (23)  vnth  Eq  (30) 


y 


rc 


+ 


■sir'-'MSi'lc  •  • 

•  firWiBlSc  ' 

crSsil'wic  • 

■  >r5  iS'^'ic' 

4 

_k=J 

5 

k^l 


5 

4  '  ~  T 

j£l  ’irji^jk  "ii^v/ic' 


■51 

<^55 


=  C  55  +  D 

where 

1  i  c  i  5 

1  s  i  s  5 

1  i  r  i  5 


(31) 


(32) 


Thus ,  Eq  (32)  is  a  ir.atrix  equation  consisting  of  250  equations  in  terms 
of  5  variables. 

Hi  n  iiiii  za  t  ion  of  Eq  (32)  accomplishes  the  goal  of  rii  ninii  zing  the 
coupling  effect.  A  cost  function  of  one  half  tin?  squarc'd  error  was 
chosen  for  this  minimization.  Thus 


J(r5)  =  ^  (C  1;^  +  0)^'  (C  ;  g  +  D) 

Sotting  the  first  derivative  of  J(,5)  equal  to  zero  and  solving  for  15 
produced 

15  -  -l.c’''cr‘  C^o  (33) 

Thus,  choosing  the  fifth  row  of  .■  (hence,  L)  in  accordance  with  Eq  (33) 
will  minimize  the  coupling  effect. 


For  the  model  matching  application,  21  sets  of  eigenfunctions  wore 


examined  in  choosing  R  (Ref  3)  with  the  best  one  being 

R  =  [162.1  25.54  -49.12  -3.17  7.26] 

This  R  assigned  the  eigenvalues  of  -jl  to  :•■]  =  -2,  X2  =  -4,  A3  =  -15, 

A4  =  -17,  and  .'.5  =  -19  with  the  eigenfunction  corresponding  to  >.3  =  -15 
being  dominant  over  the  other  ones  by  one  order  of  magnitude.  Equation 
(33)  was  employed  resulting  in 

^5''’  =  [7.89  -13.24  13.98  -4.62  -15.38] 

which  in  turn  produced 


and 


L  =  [163.41  23.33  -46.79  -3.94  2.36] 

0.99  0  0  0  0 

0  1.09  0  0  0 

K  =  0  0  1 .05  0  0 

000  0.80  0 
0  0  0  0  3.07 


In  this  application,  Gc  was  a  scalar  and  determined  to  be  3.33.  The 
above  matrices  were  substituted  into  E<)  (22)  to  produce  the  system 
state  eiiuation  given  in  Fig  24.  Time  responses  were  produi cd  with 
both  thi.  RPV  and  the  F-15  responses  on  the  same  plot  (rigs  2i)  fl). 
Thf^se  curves  sho.vcd  reasonable  agr’c.-ment ,  but  m.t  to  tiiu  desirru 
accuracy.  Ilcv.*'v.?r,  with  refincisont,  this  method  would  lU-oduce  act 
dU()l  i  ca  t  i on . 
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-0.0b9  0.33  -0.29  0  0  -0.046  0.25  -0.24 


VI.  ConcUidinq  P.e:!.’;rk^ 

COilC  1  us  i  0]iS 

An  analysis  of  the  basic  RPV  provided  the  following  conclusions: 

1.  Handling  qualities  vjere  unsatisfied  in  the  phugeid,  dutch 
roll,  and  spiral  modes. 

2.  The  thrust  vectoring  unit  had  niinimal  effect  cn  the  longi¬ 
tudinal  characteristics  and  no  effect  on  the  1 ateral -d i rectional 
characteri sties . 

3.  The  short  period  damping  ratio  and  natural  rrcgucncy,  phugoid 
damping  ratio,  and  dutch  roll  natural  frequency  varied  directly  with 
the  xcg  location.  All  other  characteristics  retrained  essentially 
constant. 

The  following  conclusions  are  based  on  the  fh\jiit  control  system 
design: 

1.  Integral  of  pitch  rate  feedback  to  the  elevator  reduced 
phugoid  oscillation. 

2.  Yaw  rate  feedback  to  the  rudder  increased  the  dutch  roll 
dampirtg  ratio. 

3.  Yaw  rate  feedback  t':'  the  ailerons  stai  ili/od  tno  spit'al  mode. 

The  model  tiiatching  application  iiiiplied  the  tollcwiirj  conclusions: 

1.  The  Gui  1 1  Ciiii n-Truxal  design  tcciinitiue  I'.rnduci'd  cn:;;p  ?n^a tors 
that  would  duplicate  only  one  state  at  a  time.  .“'Iso,  t;heso  ccir.prnsators 
wore  impractical  to  implement. 

2.  The  Fntirc  E igens ti'uc t are  Assignment  tochnique  was  not  well 


wilted  to  si ri'.i e- input  systems. 


desired  du;;]  ica tion;  ho'.vcvrr,  i:  .!-  ;  ".i  'i  d-  'jri.".-  of  accuracy. 

R  G  c  0  rs  ’on  dot  ions 

Bdsed  on  observations  iiiude  trro'rjhout  tiiis  study,  the  func.ving 
recoii'jr.ondations  are  posed; 

1.  Digital  Datcoin  is  a  very  -jeneralizcd  aleorithis  so  that  many 
different  aircraft  configurat'ions  can  utilize  ic  for  the  generation 
of  stability  derivatives.  However,  for  advanced  research  on  this  RPV 
a  more  comprehensive  algorithm  (such  as  FLEXSTA'l)  might  prove  bene¬ 
ficial  for  the  model  construction. 

?..  The  cost  function  minimized  in  the  Extension  of  Observer 
Theory  was  one  half  the  scjuared  error.  Alternate  cost  functions 
niig:  t  prcvivie  better  effectiveness  in  minimizing  the  coupling  in  the 
ert'ij'’'  equation. 

3.  Tlie  Extension  of  Observer  Theory  assumed  full  state  feedback 
Par'-.ial  state  feedback  is  a  more  realistic  assumption  vdvich  could  be 
imposed  on  this  design  technique. 
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Appendix  A 
Hand  Calculations 

Several  stability  derivatives  and  all  the  control  derivatives  were 
not  computed  by  Digital  Oatcom.  In  addition,  the  moments  of  inertia 
needed  to  be  generated.  Hence,  the  purpose  of  this  appendix  is  to 
derive  the  above  mentioned  terms.  All  notation  utilized  in  this 
appendix  is  defined  in  the  references  cited. 

Stability  Derivatives 

Seven  stability  derivatives  remained  undetermined:  Cxy.  Cx  ,  Cx^^, 

d 

Cxq.  Czjj,  Cm^J,  and  Cy^.  The  following  relationships  were  utilized  in 
their  generation  (Ref  2:19,112): 

1.  Cxu  =  -2Co  -  Uo 

The  Uq  term  is  the  compressibility  effects  and  can  be  ignored. 

Thus,  for  the  TVU  attached  cases 

Cxu  =  -2Co 
=  -2(0.05) 

=  -0.10 

and,  for  the  TVU  unattached  cases 

Cxu  =  -2(0.04) 

=  -0.08 


3. 


1^.0. 

9ot 
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Averaging  the  data  produced 


(TVU  attached) 
(TVU  unattached) 

Cx  =  0.59  -  0.37 

=  0.22 

and,  for  the  TVU  unattached  cases 

Cx^  =  0.59  -  0.35 
=  0.24 


=  0.37 

da 

1^^  =  0.35 

da 

Thus,  for  the  TVU  attached  cases 


4.  Cxq  =  0 

5.  =  -2Cl  -  Uo  ^ 

lCl 

The  Uo  term  also  represents  compressibility  effects  and  v;as  ignored. 

JU 

Thus 


Czu  =  “2Cl 
=  -2  (0.59) 
=  -1.18 


6.  Cmu  =  0 

7.  Cy^  =  0 

Control  Derivatives 

Control  derivatives  are  nondimensional  derivatives  of  force  and 
moment  coefficients  with  respect  to  the  control  surfaces.  In  all,  nine 


83 


control  derivatives  were  needed  for  the  equations  of  motion: 

^^Se’  ^’’’fie’  ^yfia’  ^"^ea*  ^Ydr’  ^*6r’  following 

empirical  relationships  were  used  to  compute  these  derivatives 
(Ref  15:  Chapters  10,  11,  and  12): 


^  •  ^Xfie  =  0 


2.  C 


-  Pi 

^66  “'^'-^f  S 


=  -(3.20)(6.38) 
30.98 

=  -0.66 


3. 


^">66 


M 

c 


Hence 


■(0.66)(11.72  - 

-  xcq) 

2.271 

,  %  mac 

10 

-2.04 

25 

-1.94 

40 

-1.84 

'^y6a  =  0 
5*  =  iagl 

=  (0.35)(0.38) 

=  0.13 

6.  Cn^a  =•-  KCl 

=  (-0.20)(0.59)(0.13) 
=  -0.015 
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Hence 


(«>5)C 


=  (5.42)(1J4)(-0.65)(1.00)CI.00)(6.97) 

30.98 


=  0.90 


8.  C 


^6r  “  ^■y6r 


Zy  cos  n  -  sin  X 


_  «  1.01  sin  0  -  (11.07  -  xcq)  sin  0 

■  13.71 

=  0.067 

„  ^  cos  a  +  Zy  sin  cx 

9-  Cn^^  -  -  - 5 - 

=  [(11.07  -  xcg)  cos  0  +  1,01  sin  0] 


xcq,  %  mac 

Cngr 

10 

-0.42 

25 

-0.40 

40 

-0.38 

Moments  of  Inertia 

The  method  utilized  in  this  computation  consisted  of  separating 
the  RPV  into  major  components  and  calculating  their  moments  of  inertia. 
Two  transformations  and  ?  summation  determined  the  overall  moments  of 
inertia  (Ref  8:  Chapter  8).  For  accuracy  purposes,  this  derivation 
was  completed  in  the  units  Ib-ir  .  .  Only  the  final  results  was 
converted  to  slug-ft^. 

Initially,  a  weight  and  balance  study  of  the  RPV's  major  components 
had  to  be  completed,  the  results  of  which  are  in  Table  XVI.  Thus, 
the  Datcom  method  yields 
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Table  XVI.  Weight  and  Center  of  Gravity 
Locations  of  the  Major  RPV  Components 


W.  lb. 

X,  in. 

Body  under  wing 

43.65 

60.74 

Boom,  left 

8.70 

90.00 

Boom,  right 

8.70 

90.00 

Engine 

14.35 

82.00 

Fuel 

30.00 

39.00 

Horizontal  Tail 

4.00 

114.45 

Instruments 

10.00 

25.50 

Lead:  Case  1 

55.83 

18.00 

Case  2 

45.99 

18.00 

Case  3 

25.20 

18.00 

Case  4 

32.28 

18.00 

Case  5 

24.18 

18.00 

Case  6 

7.07 

18.00 

Nose 

5.00 

16.00 

TVU 

20.00 

106.30 

Vertical  Tail ,  left 

5.80 

138.00 

Vertical  Tail,  right 

5.80 

138.00 

Wing 

34.00 

64.25 
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1.  Wing 


Thus 


Thus 


I  =  ^  C-Ca^  +  Cb^  +  Cc^  Cb  +  Cc  Cb^  +  Cc^] 

=  [-(5.03)3  +  (29.03)3  +  (30.25)3(29.03) 

+  (30.25)(29.03)2  +  (30.25)3] 

=  10871.69  lb-in. 2 

-  Kn  I 

^oy  '^0  I  -  Ww 
- 

=  0.70  10871.69  - 

=  1492.17  lb-in. 2 

T-  _  Wwbw2  Ki  rCr+3Ct' 
iOx  2T  L  Cr  +  Ct_ 

_  (34.00)(164.5)2  (0.86)  r30.25  +  3(24) 
24  L  24  +  30.25  _ 

=  62138.60  lb-in. 2 

%  =  lOy  +  l0;( 

=  1492.17  +  62138.60 

=  63630.77  lb-in. 2 

2.  Horizontal  Tail 

I  =  tI  ^  Cb  +  Cc  Cb2  +  Cc^] 

=  [4(15.5)3] 

=  320.25  1b-in.^ 


=  61.68  lb-in. 2 
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Cr  3Ct 
Cr  +  Ct_ 

=  (4)(60)^(0.92)  [15.5  -i-  3(15.5)1 

24  L  15.5  +  15.5  J 

=  1104.4  lb-in. 2 


24 


lOz  =  lOx  +  lOy 

=  1104.00  +  61.68 
=  1165.68  lb-in. 2 


3.  Vertical  Tail 

_  Wy  bvi2  j<5 


lOx  Y§ 


]  +  2^_^tAL  p 

rCry  +  Ctv)^ 

=  (5.8)(26)2(1.22)  f]  +  2(25.85)(15.6)  " 
18  [  (25.85  +  15.6)2_ 

=  390.49  lb-in. 2 


I  =  if  (-Ca3  +  Cb3  +  Cc2  Cb  +  Cc  Cb2  +  Cc^) 

=  “il— -  [-(19.606)3  +  (25.850)3  +  (35.206)2  (25.850) 
=  2541.91  lb-in. 2 


Thus 


lOz 


KO 


(Wy  X)2 

Wy 


0.771 


2541.91 


(113.56)2 

5.8 


245.55  lb-in. 2 


^Oy  “  lOx  ^Oz 


=  390.49  +  245.55 
=  636.04  lb-in. 2 
4 .  Body 

Initially,  the  body  was  separated  into  four  different  components  which 
needed  to  be  consolidated.  This  was  done  with  the  following  results; 
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Weight,  lb 
xcg,  in. 
ycg,  in. 
zcg,  in. 


Cases 

Cases 

1-3 

4-6 

86.05 

66.05 

74.65 

65.06 

0 

0 

-0.79 

-2.54 

Thus,  for  cases  1-3 


Tn  -  Wf_SsJ<2 

37.68 


2slq  d 


_  (86.05)(7126.11)(0.88) 
37.68 

=  74919.94  lb-in. 2 
2 


3(27.4)  .  135_ 

20351^  27.4 


In  = 

^Ox  4 


'  Ss 

_1T  tB 


{86.05)(1 .04) 
4 


7126. Il1  ^ 

nrsTT. 


=  6316.29  lb-in. 2 
lOz  =  lOy 

=  74919.94  lb-in. 2 


For  cases  4-6 


_  {66.05)(4160.45)(0.42) 
3T.68 


lO 


=  17964.73  lb-in. 2 
(66.05)(0.78) 

X  4 

=  3442.85  lb-in. 2 


3(14^47)  ,  _8: 

"2 rail  14, 


4160.45 

.rl^ 


IO2  ^Oy 

-  17964.73  lb-in. 2 


47 
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5.  Engine 


Assuming 


Assuming 


lOy  =  0.061  [|  Wpde^  +  Weie^  +  (Wp  -  We))ip2] 

=  0.061  [|  (14.35)(26)2  +  (14.35)(8)2] 

=  499.82  lb-in. 2 
iQx  =  0.083  Wpde^ 

=  0.083  (14.35)(26)2 
=  805.15  lb-in. 2 
lOz  =  lOy 

=  499.82  lb-in. 2 

6.  Fuel 

a  24  X  12  X  4  in.  tank 
lOy  =  y|-  (j-x^  +  iz^) 

=  ^  (122  +  42) 

=  400  lb-in. ^ 
lOx  =  y|-  (o-y^  + 

=  ^  (24^  +  122) 

=  1800  lb-in. 2 

7.  Instruments 

a  3  X  4  X  6  in.  pack 

lOy  =  ^  <lz^) 

lOy  =■-  ]|  [(3)2  +  (6)2] 

=  37.50  lb-in. 2 
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r 

lOx  =  v^-  V/  + 

=  [(4)2  +  (6)2] 

=  43.33  lb-in. 2 
102  =  [«x2  +  «y2] 

-  If-  [(3)2  +  (4)2] 

=  20.83  lb-in. 2 
8.  Lead 


Assuming  a  cubic  shape,  the  volume  of  lead  needed  and  the  length  of  a 
side  was 


Case 

W,  lb. 

Vol ,  in.^ 

t,  in . 

1 

55.83 

136.17 

5.14 

2 

45.99 

112.17 

4.82 

3 

25.20 

61.46 

3.95 

4 

32.28 

78.73 

4.28 

5 

24.18 

58.98 

3.89 

6 

7.07 

17.24 

2.58 

Since  a  cubic  shape  was 

assumed,  all 

the  moments  of 

inertia  will  be 

the  same. 

.0  =  "L  22 

Hence 

Case 

lo,  lb-in. 2 

1 

245.83 

2 

178.08 

3 

65.53 

4 

98.55 

5 

60.98 

6 

7.84 
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Before  the  sumniation  of  individual  inoir.onts  or  inertia  could  be 
done,  each  rtio::ient  of  inertia  needed  to  be  transformed  into  the  RPV's 
axis  system.  The  first  transformation  was  to  a  set  of  parallel  axes 
located  at  the  RPV’s  nose.  Hence 

Ix  =  S[W(y^  z2)  +  lOxl 
ly  =  ^[W(x2  +  z2)  +  lOy] 

I2  =  5:[W(x2  +  y2)  +  lOz] 

Ixz  =  ELWxz  +  lOxzl 

The  second  transformation  was  to  the  aircraft's  center  of  gravity  via 

I0x|cg  =  Ix  -  W  (7^  +  zh 

lOyicg  =  ly  -  W  (x^  +  z^) 

lOzlcg  =  Iz  -  W  (x^  +  y2) 

^Oxzlcg  -  W  xz 

Substituting  the  appropriate  values  from  Tables  XVI  and  XVII  into  the 
above  two  transformations  provided  the  RPV's  moments  of  inertia 
utilized  in  this  study  (Table  XVIII). 
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Table  X'/II.  Moments  of  Inertia  of  the 
RPV's  Major  Components 


lOx.  lb  in. 2 

lOy,  lb  in. 2 

IQz.  lb  in. 2 

l0xz>  lb  in. 2 

Engine 

805.15 

499.82 

499.82 

0 

Fuel 

400.00 

1480.00 

1800.00 

0 

Horizontal  Tail 

1104.00 

61.68 

1165.68 

0 

Instruments 

43.33 

37.50 

20.83 

0 

Vertical  Tail,  left 

390.49 

636.04 

245.55 

0 

Vertical  Tail,  right 

390.49 

636.04 

245.55 

0 

Wing 

62138.60 

1492.17 

63630.77 

0 

Body:  Case  1-3 

5316.29 

74919.40 

74919.40 

0 

Case  4-6 

3442.85 

17964.73 

17964.73 

0 

Lead:  Case  1 

245.83 

245.83 

245.83 

Case  2 

178.08 

178.08 

178.08 

0 

Case  3 

65.53 

65.53 

65.53 

0 

Case  4 

98.55 

98.55 

98.55 

0 

Case  5 

60.98 

60.98 

60.98 

0 

Case  6 

7.84 

7.84 

7.84 

0 
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Appendix  B 

Aircraft  Equations  of  Motion 

The  aircraft  equations  of  motion  are  a  set  of  nonlinear,  coupled 
force  and  moment  equations  defined  in  inertial  space.  The  purpose 
of  this  appendix  is  to  derive  these  equations  in  their  complete  form 
and  to  simplify  them  into  a  more  useful  form. 

Definitions 

Initially,  two  fundamental  definitions  needed  to  be  made,  inertial 
velocity  and  inertial  angular  velocity.  These  two  quantities  were 
defined  in  body  fixed  coordinates  as 


V  =  U  Xb  +  V  Yb  +  W  Zb  (34) 

=  P  Xb  +  Q  Yb  +  R  Zb 


Applying  Newton's  second  law  in  component  form  yielded  the  three  force 
equations. 
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Moment  Equations 

The  vector  moment  equation  was  given  by  (Ref  11:286,  418). 


™  >  K  (H) 

=  ^  (H)  + 


X  H 


where 


1 

Ix  0 

"Ixz] 

P 

H  =  [Xb  Yb  Zb] 

0  ly 

0 

j 

Q 

-Ixz  0 

Iz 

R 

Expanding  Eq  (37),  substituting  into  Eq  (36),  and  simplifying 


(36) 

(37) 


eM  =  [PIx  ^  (-R-PQ)Ixz  +  QRdz  -  Iy)]Xb 
+  [PRdx  -  Iz)  +  (P2-f  )Ixz  +  Q  Iy]Vb 
+  [PQdy  -  Ix)  +  (QR-P)Ixz  +  R  Ullb 

Recall 

eM  =  L  Xb  +  M  Yb  +  N  Zb 


Thus,  the  three  moment  equations  were 


el  =  PIx  +  (-R  -  PQ)Ixz  +  QRdz  -  ly) 
EM  =  PRdx  -  Iz)  +  {P2  -  R2)Ixz  +_Q  ly 
eN  =  PQdy  -  Ix)  +  (QR  -  P)Ixz  +  R  Iz 


Assumptions 


(38) 


Equations  (35)  and  (38)  are  the  aircraft  equations  of  motion  in 


their  complete  form.  Several  assumptions  were  made  in  their  derivation 


and  are  listed  below: 


The  mass  of  the  aircraft  remains  constant. 


-  The  aircraft  is  a  rigid  body. 

-  The  earth  is  an  inertial  reference  frame. 

-  The  Xb  and  Zb  axes  form  a  plane  of  symmetry. 
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Decoupling  the  Equations  of  Motion 

Due  to  the  coupling  terms,  the  complete  form  of  the  aircraft 
equations  of  motion  was  very  difficult  to  use.  By  assuming  perturba¬ 
tions  about  straight  and  level  flight,  these  equations  decoupled  into 
a  longitudinal  and  lateral -directional  set  of  equations. 

The  longitudinal  equations  had  only  one  control  input,  a  deflection 
of  which  would  not  cause  any  P,  R,  or  V.  Hence,  Eqs  (35)  and  (38) 
simplified  to 

eFx  =  m  (U  +  QW) 

eFz  =  m  (W  -  QU)  (39) 

m  -  Q  ly 

Similarly,  the  lateral -directional  equations  had  two  control 
inputs,  deflections  of  either  would  not  cause  any  Q.  Hence,  Eqs  (35) 
and  (38)  simplified  to 


iFy 

=  m  (V  +  RU  -  PW) 

EL 

=  PIx  “  RIxz 

(40) 

EN 

=  Pixz  +  RIz 

Longitudinal  Equations 

Simplification  of  Eq  (39)  v^as  obtained  by  assuming  stability  axes 
and  small  perturbations  about  straight  and  level  flight.  These  as¬ 
sumptions  constituted  a  linearization  about  an  arbitrary  equilibrium 
point.  Thus 

U  =  Uq  <•  u 

U  =  u 

Q  =  q 

Q  =  q 

W  =  w 

W  =  w 


9 


r 


Substituting  into  Eq  (39)  and  ignoring  second  order  terms  oroduced 
the  linearized  perturbation  equations 


eFx  =  mu 

iFz  =  m{w  -q  Uq) 
eM  =  q  ly 


(41) 


The  longitudinal  forces  and  moments  were  of  three  types:  aerodynamic, 
gravitational,  and  thrust  forces.  Assuming  a  constant  altitude,  the 
aerodynamic  forces  consisted  of  lift  and  drag  which  were  a  function  of 
U,  a,  a,  q,  and  6e-  The  gravitational  forces  were  a  function  of  0. 
Hence,  assuming  that  the  thrust  line  was  along  Xb  and  through  the 
center  of  gravity,  the  longitudinal  forces  and  moments  could  be  ex¬ 
pressed  in  functional  form  as 


5^Fx  =  Fx  (U,  a,  a,  e,  q,  6e,  6j) 

tFz  =  Fz  (U,  a,  d,  e,  q,  6e)  (42) 

liM  =  M  (U,  a,  d,  q,  <Se) 


Expanding  the  x-force  equation  of  Eq  (42)  in  a  first  order  Taylor  series 
expansion 


u  +  a  +  ^4^  a 


jU 


dO 


3a 

3q 


3a 

♦  Q  50  t 


3<‘-X  'I 


(43) 


But  ^  f.T  =  0  because  thrust  (i.e.  throttle  setting)  was  constant. 

d >x  ' 

Also 


Thus 


gravity 


-mg  sin  e 


^  -mg  cos  0 
=  -mg 
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Substituting  the  above  and  Eq  (43)  into  Eq  (41) 


mu 


.  SFy  .  3Fx  Jt- 
3q  ^  356 


(44) 


Substituting  'l  r  ^  and  the  longitudinal  stability  derivatives 
(Table  XIX)  into  Eq  (44) 


mU, 


if-  'f  -  Cxu  'u  -  ^  Cx.  i  -  Cx„  a 


-Cxq  q  +  S  9  =  Cxgg  6e 


(45) 


2Uo"^q^"  Sq 

This  same  procedure  was  used  for  the  z-force  and  pitching  moment 
equation.  Note  the  substituteion  of  W  e  UQa  was  made  in  the  z-force 
equation.  Hence  the  remainder  of  Eq  (41)  became 


-Czu  'u 


mUo  _c_ 
Sq  2U0 


Cl. 


a  +  a 


_  niMn 

2Uo  "“-q  Sq 


Cl 


q  =  Cz^g  6e 


(46) 


c  I 

■Cmu  "  pTiTT  Cm.  d  -  Cm^  a  +  q 


2Uo 


Sqc 

2Uq  ^^’’q 

Equations  (45)  -  (47)  needed  to  be  rearranged  into  a  more  useful 
form,  i.e.,  a  state  equation.  Define 

Cl  =  ^  Cl 

Sq  2Uo  ‘‘d 
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Table  XIX.  Definitions  of  the  Longitudinal  Stability  Derivatives 
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Solving  Eq  (46)  for  a 


•  =  ^Zli  I,.  .^Lg  +  ] 

“  irr  "  “CT  Cl 


mUo  _c_  p 
Sq  ■  2U0  Cq 


q  +  «e 


Substituting  Eq  (48)  into  Eq  (45)  and  solving  for  'u 
I-  _  f  So"  p  ,  S^  Cxa  Czu"]  , 

[  iiife  ^  2iu52  cr^J 

Cy  +  -^-3.0^  Cxg  Clgl 


mUo 


-_3. 

UO 


,  Smr  Cx.  mU£[  cCLq  .  P^  , 

L^V  -Ct  “Sq  ■  ^[l■S  "  2^052 

Cx^Cz6e  4.  %  p„  1  , 

2^  Cl  iUn 


2mUo^  Cl  ■  mUo 
Substituting  Eq  (48)  into  Eq  (47)  and  solving  for  q 


q  = 


c„„  t  sod 

2UoIy 

Sq^  r  Sgo^^  Cm, -Cl 

-  2U0ly 


(48) 


(49) 


-_2 

Stic. 


2UoIy  Cl 


niUo _ ^  p. 

■  2Uo  ^Cq 


S_qc  Cmjg  +  Sqc 


-r2 


Cm.' C  z.ij^ 


L  ‘y 


2UoIy 


Cl 


2uoiy  ^"^q 


50 


(50) 


Choosing 


X]  =  'u 
X2  =  a 

X3  =  e 

H  =  q 
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Eqs  (48)  -  (50)  formed  the  desired  state  equation  as 


where 


A  =  A  X  +  B  Secmd 


Cl 

All 

Ai2 

Ai3 

A14 

A2I 

A22 

A23 

A24 

A3I 

A32 

A33 

A  34 
A41 


mUn  , 
Sq 


c 

2Uo 


Cl. 


Sg_ 

mUo 


Cxu 


+ 


2niL^ 


CxaCzu 

Cl 


Sq  ^  S-qc*  Cx^CL„ 
mUo  -  TmU^  Cl 


UO 


Sqr 

2mUo 

Czu 

Cl 


Cx^  mUp 
Cf  Sq 


2Ud 


Clq 


0 

0 


1 

S^  P  Sq^  Cm^Czu 

“ly  2U0ly  "  Cl 


(51) 
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ly 


2U0ly  Cl 


A43  =  0 

"  lUoTy 


C|T1‘ 

Ct 

Cl 


mUo 


ni 


mUo 


C  + 

2mUo^ 


C  CLq 

■  -w 

CxaCz^e 

Cl 


+  c 


mq 


B21  = 


'Z6e 


831 

B41 


=  0 


ly  "'"5e  2U0ly 


Cl 


Lateral -Pi rectional  Equations 

As  in  the  longitudinal  equations,  a  linearization  about  an 
arbitrary  equilibrium  point  was  completed.  Thus 


U  =  Uq  +  u 

V  =  V 

V  =  » 


W  =  w 


P  =  p 
P  =  p 
R  =  r 
R  =  f 


Substituting  into  Eq  (40)  and  ignoring  second  order  terms  produrod  the 
linearized  perturbation  equations 


tFy  =  m(v  +  llQr) 

::L  =  p  Ix  -  f  Ixz 
tN  =  -p  Ixz  -  r  Iz 


(52) 
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The  latcral-df rectional  forces  and  moments  consisted  only  of  aero¬ 
dynamic  and  gravitational  forces.  As  before,  the  aerodyiamic  forces 
were  lift  and  drag;  however,  in  this  case  they  were  a  function  of  3, 
r,  p,  5a,  and  6r.  The  gravitational  forces  wore  a  function  of  only  f. 
Hence,  in  functional  form 


5:Fy  =  Fy  (3,  r,  p,  6a,  6r,(j)) 

tL  =  L  {3,  r,  p,  6a,  5r)  (53) 

iN  =  N  (b,  r,  p,  6a,  6r) 


Expanding  the  y-force  equation  of  Eq  (53)  in  a  first  order  Taylor  series 
expansion 


Recal 1 

Therefore 


Also 


r  +  p 
3p 


+  Ifl 

36a 


6a 


+  iEz 

36r 


sr  + 


sFy 

3<i. 


4- 


''J'lgravity  "  ”3  ♦ 

^  =  mg  cos  ♦ 

=  mg 


V  r:  Uq 
v  •  UqB 

Substituting  the  above  and  Eq  (54)  into  Eq  (52) 


(54) 


mUo  (’>  +  r)  =  3  +  r  +  p 

H'.  jr  op 


+  6a  +  6r  +  mq  ■}: 

.^ca  o6r 


(55) 
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Substituting  the  idteral -directional  stability  derivatives  (Table  XX) 
into  Eq  (55) 


P  mUo 

-Cy3  ^  e  + 


r.i'-'o 


S^q  ■  2Uo 


r  - 


mg 


6r 


(56) 


This  same  procedure  v;as  used  for  the  rolling  and  yawing  moment 
equations  with  the  following  results 


■^"6  ^  ■  2Uo  ■  2Uo  ^“^P  ^  ^  SqT)  ^ 


=  5a  +  6r 


S  -  2^jQ  Cnr  r  +  5^  r  -  ^Uq  ^"P  P  'I^o  P 


-  Cn^3  fa  +Cn^^  6r 


(57) 


(58) 


As  before,  Eqs  (56)  -  (58)  needed  to  be  put  in  state  equation 
form.  Solving  Eq  (56)  for  3 


e  =  4  cyf. 


2rnUo‘^ 


1 


r  +  Ff^ 

Uo 


+  Cv„  P  +  :::rr^  fa  +  Cutv. 


2mUo2  P  mUo  ^oa  mUo  ^^r 
Solving  Eq  (57)  for  p 


(59) 


^  c 

.  +  ,.Ssbl 

Co  r  + 

ix  z,  ^ 

Ix 

2UoIx 

^  \  ' 

Ix 

SjLb-  r 

.  p  + 

P  0  f 

+  S|i 

2UoIx  ' 

P  P  Ix 

^■«a  “ 

Ix 

(60) 


Substituting  Eq  (60)  into  Eq  (58)  and  solving  for  r 
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Table  XX.  Definitions  of  the  Lateral-Directional  Stal)ility  Derivatives 


1 

cl 

Sqb  p  , 
^  Iz  "S  ^ 

Ixz  ^b  p 

1 

Cn 

.  Ixz 

Sgb2 

C2‘ 

L  2UoIz 

Iz 

■2UbIx 

1 

Cnn 

.  Ixz 

Sqb^ 

C2 

2UoIz  P 

Iz 

2U0lx 

1 

Cn 

+  ^ 

Sqb 

C2 

Iz  "6a 

Iz 

^Ix 

1 

Cn 

+ 

Sqb  , 

C2 

Iz 

Iz 

Ix 

■] 


6d 


(61) 


where 


C2  =  1  - 


Ixz^ 
Ix  Iz 


Substituting  Eq  (61)  into  Eq  (60)  and  solving  for  p 


P  = 


Ixz' 


Wiz  ^"3  ^ cnYTY 


51^2  Ixz 


Ixz  Tn  +  Sqb^  P  ■,  ^ 
_2U^  Ix  Iz  2UoIx  L 


C2 


Ix 

Iz 

Ixz 

Ix 

Iz 

Ixz 

Ix 

T7 

Ixz 

Ix 

T7 

Cnr  + 


Ctr 


5'qb2 

2UoIx 
2 


1  + 


Ixz‘ 


C2  Ix  Iz 
Ixz^ 

C2  Ix  Iz 


1  + 


Ixz^ 


'6r 


6r 


C21xlz 
2 


1  + 

•  ~  ro  1 


C21xlz 


6a 


6r 


(62) 
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Choosing 


=  a 


Xo  = 


X4  =  p 

Eqs  (59),  (61),  and  (62)  formed  the  desired  state  equation  as 


where 


C2  =  1 


Ixz^ 

Ixlz 


Cj,,  -I 

*13  =  uf 


*'4  “ 


^21  " 


C2 


^  Cn  +  5gb  r„ 

I2  ^'^c<  I-’  T'' 


Ix 


Apo  =  -  ’  *"5^  ^  Ixz  sn^.  c 
^22  c.  ,  2UoIz  Iz  2UoIx  i" 


A23  =  0 


Mi 

2UoIz 


+ 


Sqb^ 

2Uorx 


(63) 


108 


A31  =  0 


A32  =  0 


A33  =  0 


A34  =  1 


A41 


_  Sgb  Ixz_  r  +  Cp 
■  C2  I7I2  lx 


1  + 

I  ^  ro  T 


1X2^ 


C21xlz 


^42  =  1^7  Cnr  +  2U^ 


1  +  Ixz^ 
C2l71z 


A43  =  0 


A 


Sqb^ 
44  "  2UoC2 


1x2 

ixlz 


Cnp  ■*■ 


Sqb^ 

2UoIx 


1  + 


Ixz?_ 

ZTUlz_ 


1  =  iDo  ^yaa 
^12  =  iil§  ^y6r 


B21 

1 

■  C2 

'sgb 

L  Iz 

Sgb 

Ix 

^22 

1 

"  C2 

[Sgb 
_  Iz 

Cnar  *  ^ 

Sqb 

■^x 

- 

B3I 

=  0 

B32  =  0 

B4I  =  ^  ^ 


6  a 


1  +  _IX2_ 
'  C21xlz 


B42 


-  Sgb  Ixz  r  +  Q. 
C2  Ixlz  Ix  " 


<sr 


1  + 


J_xzi_ 

C21xlz 
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Summary 

The  complete  form  of  the  force  and  moment  equations  can  be  found 
in  Eqs  (35)  and  (38).  The  decoupled  linearized  perturbation  equations 
can  be  found  in  Eqs  (41)  and  (52).  Lastly,  the  longitudinal  and  lateral- 
directional  state  equations  can  be  found  in  Eqs  (51)  and  (63). 


no 
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Appendix  C 
Time  Responses 

Representative  time  responses  were  generated  for  the  basic  RPV, 
the  augmented  RPV,  and  the  nonlinear  simulation  and  are  contained  in 
this  appendix.  Each  figure  contains  a  series  of  three  time  responses, 
one  for  each  of  the  systems  mentioned  above.  This  arrangement  facili¬ 
tates  easy  comparison  between  the  different  systems.  The  longitudinal 
time  responses  are  displayed  in  Figs  30  -  35.  The  lateral  time  re¬ 
sponses  are  given  in  Figs  36  -  44.  Lastly,  the  directional  time 
responses  are  contained  in  Figs  45  -  53. 


Ill 


u 

<u 


3 


20 


Nonlinear  Simulation 


I  / 
V' 


0  ^  2o"  40  '  60  "^80  ^  foo 


Time  (sec) 


Figure  30.  RPV  Forward  Velocity  due  to  a  1  sec  Pulse 
Elevator  Command  of  -5  dog 
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Figure  31.  RPV  Angle  of  Attack  due  to  a  1 
Pulse  Elevator  Command  of  -5  deg 
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Figure  33.  RPV  Pitch  Rate  due  to  a  1  sec 
Pulse  Elevator  Command  of  -5  deq 
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36.  RPV  Sides!  i;,'  A'vju'  due  to  a  1 
Pulse  Aileron  Cosjiand  of  5  deg 


Fitjurn  37.  RPV  Yav;  Anqle  duo  to  a  1  sec 
Pulse  Aileron  Coniniand  of  5  deg 
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Figure  38.  RPV  Yaw  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  39.  RPV  Bank  Angle  due  to  a  1  sec 
Pulse  Aileron  Convnand  of  5  deg 
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Figure  40.  RPV  Roll  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  42.  RPV  Aileron  Deflection  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  43,  RPV  Rudder  Deflection  Angle  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  44.  RPV  Rudder  Deflection  Rale  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 


Figure  45.  RPV  Sideslip  Angle  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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Figure  46.  RPV  Yaw  Angle  due  to  a  1  sec 
Pulse  Rudder  Comir:and  of  5  deg 
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Figure  48.  RPV  Bank  Angle  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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Figure  50.  RPV  Aileror.  Deflection  Angle  due  to  a  1 
Pulse  Rudder  Coirmand  of  5  deg 
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Figure  52.  RPV  Rudder  Deflection  Angle  due  to  a 
Pulse  Rudder  Command  of  5  deg 
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Figure  53.  RPV  Rudder  Deflection  Rate  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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flight  control  system  so  that  the  RPV  has  dynamic  characteristics  similar  to 
the  F-15. 
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